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Introduction 



The important topic of lattice renormalization and improvement, covered by P. Weisz 
in this School, is quite complicated, even for students who have considerable exper- 
tise in calculations based on continuum regularization (e.g. Pauli-Villars, dimensional 
regularization). The reader is advised to familiarize himself with the basic concepts, 
before tackling the present chapter, which is meant to be a complement to Chapter 2 
of this book ( [Weisz, 2010p . Our aim is to present three very specific topics: 

• The consequences of the loss of chiral symmetry in the Wilson lattice regulariza- 
tion of the fermionic action and its recovery in the continuum limit. The treatment 
of these arguments involves lattice Ward identities. 

• The definition and properties of mass independent renormalization schemes, which 
are suitable for a non-perturbative computation of various operator renormaliza- 
tion constants. 

• The modification of the Wilson fermion action, by the introduction of a chirally 
twisted mass term (known as twisted mass QCD and abbreviated as tmQCD), 
which results to improved (re)normalization and scaling properties for physical 
quantities of interest. 

Loss of chiral symmetry implies that the quark mass is not multiplicatively renor- 
malizable, having also an additive counter-term, proportional to the ultraviolet cutoff 
(i.e. the inverse lattice spacing). Moreover, the no-renormalization theorem of the par- 
tially conserved vector and axial currents (PCVC and PCAC) turns out to be some- 
what intricate. Wilson fermions preserve flavour (vector) symmetry, with an exactly 
conserved current, which is a point-split discretization of the corresponding contin- 
uum operator. The vector local current Vp,{x), familiar to us from continuum QCD, is 
to be normalized by a scale independent normalization factor Zv{go) 7^ li which is a 
function of the bare coupling go, tending to unity in the continuum limit. On the other 
hand, axial symmetry is broken by the Wilson term. Thus any discretization of the ax- 
ial current on the lattice is accompanied by its scale independent normalization, which 
tends to unity in the continuum limit. For instance, the familiar axial current A^{x) 
has a normalization ZA{go) 1- Another consequence of the loss of chiral symmetry 
is that operators belonging to the same chiral multiplet do not renormalize by the 
same renormalization factor. Although they all have identical anomalous dimensions 
(i.e. the same scale dependence), ratios of their renormalization constants arc scale 
independent quantities which become unity only in the continuum limit. These prop- 
erties are established with the aid of Ward identities. The very same Ward identities 
may also be used in practical computations for the non-perturbative determination of 
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these scale invariant current normalizations and ratios of rcnormalization constants, 
at fixed lattice spacing. 

ft is a well established fact that lattice perturbation theory is characterized by 
poor convergence. Renormalization constants calculated perturbatively are therefore 
sources of rather large systematic errors in lattice predictions and postdictions of many 
interesting physical quantities. The second topic of these lecture notes is the discus- 
sion of a family of renormalization schemes (RI/MOM, RI/SMOM) which are used for 
the non-pertubative renormalization of fermionic composite operators. The two main 
sources of systematic error (unwanted low energy effects in the infrared and discretiza- 
tion effects in the ultraviolet) necessitate the existence of a "renormalization window" , 
in which these errors are controlled. The infrared problems due to a Goldstone pole 
in some renormalization constants is discussed in considerable detail. Although our 
presentation is based on Wilson fermions, these renormalization methods arc equally 
well applied in other lattice regularizations, such as domain wall fermions. 

The last topic concerns the benevolent effects of tmQCD in the rcnormalization and 
improvement properties of physical quantities. For this modified Wilson fermion action 
the Wilson term (responsible for a hard breaking of chiral symmetry) and the twisted 
mass term (which breaks the symmetry softly) are in some sense orthogonal in chiral 
space. This results to some fermion operators having simpler renormalization patterns 
in tmQCD than in the standard Wilson lattice theory. In these circumstances, the 
computation of their matrix elements is under better control in the tmQCD framework. 
Moreover, discrete symmetries in tmQCD imply that a particular tuning of the twisted 
mass action results to 0{a) improvement of physical quantities, without the need of 
Symanzik counter-terms. This is the so-called "automatic improvement" of Wilson 
fermions. 

In these lecture notes we have strived to give explicit derivations of the most 
important results. Having opted for a detailed treatment of the theoretical issues, we 
have not addressed the many interesting numerical results, which have appeared in the 
literature over the years. Moreover, since these are meant to be paedagogical lectures 
addressing this advanced subject at an introductory level, we have only included the 
essential references which would facilitate the student's further study. This inevitably 
introduces some bias, for which we apologize. 
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Basics 



Although we understand that the student is famihar with the basics of quark mass 
and composite operator renormalization, we recapitulate them here for completeness 
and in order to fix our notation. We also collect several useful definitions; although 
this is somewhat tedious, it is important to spell out the not-so-standard notation 
right from the beginning. 

Concerning notation, we have preferred economy to mathematical rigour. Since 
the various bare quantities discussed below are defined in the lattice regularization, 
integrals like say, J d^xid'^X2 of eq. (|2.28p . are really sums (a® ^^), which run over 
all lattice sites, labelled by xi and X2 etc. The occasional use of integrals instead of 
sums, partial derivatives instead of finite differences, and Dirac functions instead of 
Kronecker symbols, simplifies notation, hopefully avoiding any confusion. Moreover, 
a space-time function and its Fourier-transform will be indicated by the same symbol 
with a different argument (e.g. f{x) and f{p))\ again mathematical rigour is being 
sacrificed in favour of notational economy. 

2.1 The Wilson lattice action and its symmetries 

We opt for the lattice regularization scheme, proposed by Wilson, which consists of a 



gluonic actioiij([Wilson, 19741 



5o V L ^ 



(2.1) 



and a fermionic one ( [Wilson, 1977[ ), 

X 

^ 2a 



Mo 



■ip{x) 



i^{x){l ~ 7^)[/,,(x)V'(x + afi) + ^{x + aA)(l + j^)Ul{x)^{x) 



V'(x) 



(2.2) 



In standard notation, Ufj,{x) = cxp[igoaGfj,{x)] is the lattice gauge link, depending on 
the gauge (gluon) field Gi_i{x), go is the bare coupling constant, a the lattice spacing 

Clearly renormalization schemes, constants etc. depend on the lattice actions. However, the choice 
of a specific gluonic action, such as that of eq. 1 12. Il l, is of no particular consequence to the topics 
discussed in the present lectures. 
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(i.e. our UV cutoff) and Up the Wilson oriented plaquette. Two discretizations of the 
covariant derivative are used in the above definition: 

aV^.'tpix) ^Uf,{x)-)jj{x + afi) -ipix) (2.3) 
aVf^ ip{x) = ipix) - U'^{x - afi)ip{x - afi) . (2.4) 

Later on we will also make use of the backward covariant derivative, defined as: 



aip{x) I?p= i){x + afL)Uj^{x) - ijj{x) 



(2.5) 



The quark field il'{x) is a vector in flavour space. Its components are denoted by 
(/ = 1, 2, 3, . . . , A^f)- The diagonal bare mass matrix is denoted by Mq and its 
elements by mo/ (/ = 1, . . . ,Np). Two continuous symmetries of particular interest 
are the flavour and chiral symmetries. The flavour group SU{Nf)v consists in the 
global vector transformations of the fcrmion fleld: 



V'(x) 
ip(x) 



ip'ix) 

V7(x) 



exp 



tp{x) exp 



tpix) 
„ A° 



-I a. 



(2.6) 



with a — 1, 



, Np — 1. The group generators in the fundamental representation are 



the flavour matrices A''/2 (a = 1, . . 

Tr[A'^A''] 
A'' A''' 

~2' y 

A^ A^ 
T'Y 



,iV| - 1), satisfying: 



ab 



26 



■' 2 



abc ' 



Nf 



(2.7) 



where / represents the Np x Np unit matrix in flavour space; /"^'^ are the SU{Np) 
structure constants and d"''^ are totally symmetric. For Np = 3 the generators A° are 
the eight Gell-Mann matrices, while for Np = 2 they are the three Pauli matrices. 
Under these transformations the Wilson lattice action (|2.2p is invariant, provided all 
quark masses are degenerate (mgi = • • • = niQNp). Vector (or flavour) symmetry is 
preserved by the Wilson lattice regularization. 

Next we consider the axial transformation of the fermion field: 



ip{x) 
ip{x) 



^\x) 



exp 



%j}{x) 



ipix) exp 



i a 



A" 



(2.8) 



with a — 1, • • • , Np — 1. Even in the absence of a mass matrix {Mq = 0), the presence of 
the Wilson term {aip V^i^Vfj, V') in the action (|2.2p ensures that it is not invariant under 
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these transformations. Consequently, the chiral group SU {Np) S'[/ (A^f)^, related to 
axial and vector transformations, is not a symmetry of Wilson fermions. An interesting 
generalization of chiral symmetries is described in Appendix [71 For completeness, we 
also collect the discrete symmetries of the Wilson action in Appendix [51 

We close this subsection with a general observation. Many important quantities 
we will be using, such as the quark propagator and the various correlation functions 
Go, Ao, To (to be defined below), are gauge dependent. This implies that some gauge 
fixing procedure has been applied in their calculation. The gauge of choice is usually 
the Landau gauge, which on the lattice is understood as imposing a discrete version of 
the condition c^^G^ = on all lattice gauge fields G^. The reader is advised to consult 



a review of lattice gauge fixing for more details ( jGiusti et ai, 2001 1 . 

2.2 Quark mass renormalization 

We recapitulate the general renormalization properties of the quark masses with Wil- 



son fermions. It is convenient ( Bhattacharya et ai, 20061 to separate non-singlet and 



singlet quark mass contributions of the bare mass matrix: 

Mo = J^d^"^^"^ + ^ ■ (2.9) 

The sum runs over diagonal group generators only (e.g. d = 3 for SU (2)v and c? = 3, 8 
for SU{3)v) and m^^ is the average of the diagonal elements of Mq. Note that in 
the flavour symmetric theory, where all masses are degenerate, we have fh'^ = and 
Mq = m^^I; cf. eq. (|7.7p . In terms of the mass decomposition (|2.9p . the bare mass of 
a specific quark of flavour / is 

mo/ = m^^ -I- ^^m'^Xjf (2.10) 

with X'jj the /"^ diagonal entry of the matrix A''. 

As shown in Appendix [71 spurionic flavour symmetry implies that all components 
m° transform as a multiplet in the adjoint representation of SU{Np)v , while m'^^ is a 
singlet. Thus, as far as this symmetry is concerned, members of the adjoint multiplet 
renormalize in the same way, while the singlet m^^ could renormalize differently. In the 
same Appendix we have also seen that spurionic axial transformations mix flavour non- 
singlets m° and singlets m^^ . If that were a symmetry of the theory, then and m^^ , 
being members of the same chiral multiplet would share a common renormalization 
constant. This is not the case, however, as the Wilson term of the lattice action breaks 
chiral spurionic symmetry. Moreover, flavour symmetry also suggests that the fh'^ 
components renormalize multiplicatively, while the lack of chiral symmetry allows 
m'^^ to mix with the identity. Thus the singlet mass m**^ is also subject to additive 
renormalization. The bottom line is the following renormalization pattern: 



" a— !-0 



[m-^(,9R,)]p = hm Z„o(5o,aA^) [ m-^(go) - m„. ] . (2.11) 

We denote by Z^n, Z„fi the multiplicative renormalization constants and ^ the renor- 
malization scale; gn^n) is the renormalized gauge coupling. One must also keep in 
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mind the dependence of the bare coupling on the lattice spacing; go (a) ^ l/ln(a). 
The additive renormalization factor mcr = [C{go)/a\, with C{go) a coefficient bear- 
ing the dependence on the bare coupling, is a flavour independent, linearly divergent 
counter-term. At tree- level, mcr = — 4/a. This counter-term is a well known feature of 
Wilson fermions. As already pointed out, it is a consequence of the loss of chiral sym- 
metry by this regularization (even when the bare quark masses are switched off). The 
fact that there are two mass renormalization constants, Zm 7^ -Z^m"! albeit with the 
same scale-yit dependence, is also due to the loss of chiral symmetry. We will establish 
useful relations between them in sect [3l 

It is common practice to omit, in expressions like (|2.1ip . the continuum limit on 
the r.h.s.. The resulting equations relate quantities at fixed lattice spacing and are 
true up to discretization effects. 

Given that eq. (|2.10p relates the bare mass with flavour / to the singlet and non- 
singlet bare masses, we define the renormalized mass of this quark as the same combi- 
nation of the renormalized singlet and non-singlet masses ( [Bhattacharya et ai, 2006| : 

= Z„o[m^^(go) - mcr] + ""'^ -^/Z 

= Z„o[m""(go) - mcr] + ^m[mo/ - m^^ . (2.12) 

The second line is obtained by using eqs. ()2.11|) . The third line is obtained from the 
second, by substituting the sum in its last term from eq. (|2.10p . This result can also 
be conveniently rewritten as 

[m/]R = Zm 

Note that it is not enough to take mo/ — > mcr in order to obtain the chiral limit 
[m/]R — > 0; this is achieved only when all bare quark masses moi-, ■ ■ ■ moNp tend 
to mcr- Also note that if the regularization scheme were chirally invariant (it isn't!) 
then mcr would be absent and Zm = -^mo; in this case the quark mass would be 
multiplicatively renormalizable ([m/]R = Z^mo/), for each flavour. 

It is instructive to consider the case of degenerate masses {rh'^ = and m^^ = mo/). 
Then eq. p.l3p reduces to the more familiar 

[m/(.gR)]R = lini[ Z„,o(go,a^) m/(.go,mo/) ] , (2.14) 

a— >0 

where the subtracted quark mass is defined as 

m/(go,mo/) = mof{go) ~ m„ . (2.15) 

The chiral limit is now simply m/ — 0; i.e. mo/ mcr- In the mass degenerate case, 
or whenever flavour dependence is unimportant, the flavour index will be dropped, 
leaving us with mo, m, mR etc. 



mo/ 



mcr) 



(2.13) 
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We also point out that combining eq. p.l3p for two flavours, say / = 1, 2 we obtain 

[ mi ]r - [ m2 ]r = Zrn[moi - too2 ] 

= Zrn[mi - 7712] ; (2-16) 

i.e. bare (and subtracted) mass differences renormalize with Z,„ rather than Z,„o. As 
already anticipated, the two renormalization constants and Z„jO have the same 
scale dependence. We will show this in sect. [31 and also prove that their ratio is a 
finite function of the bare coupling go, which tens to unity in the continuum limit. 



2.3 Quark propagator 

The bare quark propagator in coordinate space, Sf{xi —X2] m) = {'ipf(xi)'>pf{x2)), has 
the Fourier Transform 



S{p;m) = d X exp{~ipx) S{x,m) , 



(2.17) 



where the flavour index / has been dropped (it will reappear wherever necessary) . The 
renormalized propagator is then given by 



5(p;gR,mR,/x) ]r 



lim [ Z^{go,afi) S{p;go,m) 

a— 5-0 



(2.18) 



1 /2 . . . . 1/2 

where Z^ is the wave function renormalization; ipn, = Z^ ip. Note that in the bare 
quantities, the bare mass mp has been traded off for the simpler (from the point 
of view of the chiral limit) subtracted mass m. For later convenience, we define two 
"projections" (i.e. traces in spin and colour space) of iS~^(p, m) and its derivative with 
respect to momentum: 



rs(p;m) = ^ Tr 



Tmip; m) 



12 



Tr 



dS ^ {p; m) 
dp, 
dS~^{p] m) 



dm 



(2.19) 
(2.20) 



These quantities have been defined so that their tree level values are equal to unity. 
From eqs (j2.14p and (|2.16p we see that T^{p) and Tm.{p) renormalize like 



[rs(p)]ri = 1™ ^^/rs(p,m) 

[ rm(p) ]r 



lim 

a-i-O 



(2.21) 
(2.22) 



It is natural to fix the quark field and quark mass renormalization Z^ and Zm by 
imposing renormalization conditions on [rE(p)]pj and [rm(p)]j^ at a given momentum 
p = H- Such a scheme will be the subject of Section 21 It is clearly not a unique choice. 
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2.4 Quark bilinear operators (composite fields) 

In these lectures we will make extensive use of local bilinear quark operators of the 
form: 

Olix) = {^{x)T^^P{x) , (2.23) 

where T stands for a generic Dirac matrix and / = 1, • • • , Np — 1 (flavour non-siglet 
case). Specific bilinear operators will be denoted according to their Lorentz group 
transformations: the scalar, pseudoscalar and tensor densities are 

S^{x) = ■ilj{x)—ij{x) , P'{x) = 'ip{x)—-f5ij{x) , 

Tf,{x) = ^(x)^7^7,^(a:) , (2.24) 
respectively, whereas the vector and axial vector currents are 

V^ix) = ^ix)^^^i,{x) , Alix) = ^(a;)^7^75^(x) . (2.25) 

In the above, implicit colour and spin indices are contracted. Besides the non-singlet 
bilinear quark operators defined in eqs. (|2.24p and (|2.25p . we will also be using the 
singlet ones 

S°ix) = ^(x)y^(x) , P°(x) = ^{x)^j54ix) , (2.26) 



with A° = Nf I ■ 

We will sometimes simplify the notation by using the symbols 5", etc. (dropping 
the flavour superscript /). These are meant to be non-singlet operators of the form 
Or = 4'i^'4'2- For example if we substitute A-^ with A+ = Ai + 1X2 in eqs. p.24p 
and ()2.25p . we end up with, say, ud, u"ffj,d etc.. The renormalization and improvement 
of four-fermion operators, although very important, is beyond the scope of the present 
lectures. 

The insertion of the operator Or = iJiT^p2 in the 2-point fermionic Green function 
gives the vertex function 

Goixi-x,X2-x) = {Mxi)Orix)Mx2)) , (2.27) 

with translational invariance explicitly taken into account in the notation. Placing, for 
simplicity, the operator at the origin a; = 0, the above correlation function becomes, 
in momentum space, the vertex function 

Go{pi,P2) = J d'^xid'^X2 exp{-ipixi) exp{ip2X2) Go{xi,X2) ■ (2.28) 
The corresponding amputated Green function is given by 

Ao(pi,P2) = S^'ipi) Go{pi,P2)S2'ip2) , (2.29) 

where the propagator subscripts denote flavour. Note that Go and Aq are rank-2 
tensors in colour and spin space; the colour and spin indices are suppressed in the 
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notation. Since it is preferable to work with scalar, rather than tensor-like correlation 
functions, wc define the projected amputated Green function Toip) as follows: 

ro(p) = ^Tr [PoAo{p,p)] . (2.30) 

The trace is over spin and colour indices. In some cases Lorentz indices need to be 
added; e.g. the axial current has vertex functions G^, etc. Thus, with the 
normalization factor 1/12, Pq is the Dirac matrix which renders the tree-level value 
of To{p) equal to unity; i.e. it projects out the nominal Dirac structure of the Green 
function Ao{p)- It is easy to work out these so-called "projectors" for each of the 
operators of eqs (PTM]) and (j^^ : 



1 „„ 1 



- ; = 4757m • (2.31) 

Repeated Lorenz indices, appearing in the definition of these projectors and the cor- 
responding Aq's are meant to be summed over. For simplicity of notation, specific 
Green functions will be denoted, for instance, as Gs,Fp, Ay etc. 

The dimension-3 bilinear operators S-^,P-f and Tj^^ (with a = f,...,Np — 1), 
are non-singlet in flavour space. As there are no other operators of equal or smaller 
dimension with the same quantum numbers, they renormalize multiplicatively; their 
rcnormalization constants arc denoted as Zs,Zp and Zt- Thus, the rcnormalizcd 
operator is formally given by 

[ Or(5R,"iR.Ai) ]r = ^^m [Zo{go,afj,)Or{go,m)]. (2.32) 

a— >0 

For the currents V/ and A(, there is a lattice version of the no-rcnormalization theorem, 
which will be the object of the following sections on Ward identities. Essentially we 
will show in Sect. [3] that these currents, once normalized by appropriate finite, scale 
independent factors Zy and Za, go over to their continuum limit up to discretization 
effects which arc proportional to positive powers of the lattice spacing. We then write 
the continuum currents as: 

[ Ov,Ai9K,mR) ]r = lim [ Zv,a{9q) Ov,A{ga,m) ] ■ (2.33) 

Combining the above wih cqs. ()2.18|) and ()2.27|) - ()2.29|) . wc find for the rcnormalizcd 
Green functions: 

[ Go{pi,P2) ]r = lim [ Z,^, Zo Go{pi,P2) ] , (2.34) 
[ Ao(pi,P2) ]r - lim [ Zo Ao{pi:P2) ] • (2.35) 
The renormahzation of To{p) is identical to that of Ao(p): 

[ To{p) ]r - lim^ [ Zo To{p) ] ■ (2.36) 



a-i-O 



The singlet pscudoscalar density P*^ renormalizes multiplicatively in a similar fashion; 
its renormahzation constant is denoted by Zpo. The singlet scalar density 5°, besides 
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its multiplicative renormalization constant Zgo , also requires a power subtraction, as 
we will discuss in section [XH 

The above definitions of renormalization constants are purely formal. By this we 
mean that some arbitrary (still unspecified), choice of renormalization scheme and 
scale is implied. No reference to a specific scheme has been made so far. However, it is 
implicit in the notation that only mass independent schemes are taken into consider- 
ation (otherwise Zm, and Zq would also depend on quark masses). In sect. 14.11 we 
will discuss a specific mass independent renormalization scheme, known as RI/MOM 
dMartinelh et al, 1995|. 



3 

Lattice Ward identities 



Ward identities are relations between Green functions; essentially they express how a 
classical continuum symmetry is realized at the quantum level. Since lattice actions 
break chiral symmetr}0 the resulting Ward identities are not trivial discretized tran- 
scriptions of the formal continuum ones. In fact they dictate the proper normalization 
of (partially) conserved currents on the lattice and relate the renormalization param- 
eters of quantities which in the continuum belong to the same chiral multiplet, such 
as the scalar and pseudoscalar densities. In this way the question of recovering chiral 
symmetry in the continuum limit (up to discretization effects generically denoted here 
as O(a)-efFects) is linked to the appropriate (re)normalization of lattice operators in 
the lattice regularization. In particular we will discuss the following topics in detail: 

• the conservation of vector symmetry SU{Np)v on the lattice with Wilson fermions; 

• the loss of chiral SU{Nf)l^SU {Nf)r symmetry on the lattice, even in the chiral 
limit, and its recovery in the continuum limit; 

• the derivation of vector and axial lattice Ward identities and the finite normal- 
ization of the vector and axial currents, resulting from them; 

• the relations between quark mass renormalization and the renormalization of the 
scalar and pseudoscalar densities, arising from lattice Ward identities. 

Let us first recall that at the classical level. Ward identities connect the four- 
divergence of vector and axial currents to the scalar and pseudoscalar densities. The 
so-called "partially conserved vector current" (PCVC) relation 



results, through the standard Noether construction, from the invariance of the QCD 
(Euclidean) fermionic action under SU{Nf)v transformations. Vector current conser- 
vation holds in the degenerate mass case (Mq (x /). The "partially conserved axial 
current" (PCAC) relation 



^This statement is usually reserved for the Wilson fermion action only; staggered fermions are 
known to display a reduced chiral symmetry, while Ginsparg- Wilson ones preserve a lattice chiral 
symmetry. Yet, in practice even with the latter actions there is some loss of chirality: for domain wall 
fermions the extension of the fifth dimension is never quite infinite. 



d^V^ix) + ^(x) — ,Mo ^J{x) = , 



(3.1) 




(3.2) 
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results from the invariancc of the QCD action under axial transformations. Axial cur- 
rent conservation holds in the chiral limit (Mq = 00 Similar, possibly more familiar, 
expressions are obtained in the specific case with A" replaced by A+ = A^ + This 
involves transformations on two flavours only, tpi and '02 (with corresponding masses 
mi and 7712). With this choice, the quantities S, P, T^jy, V^, of eqs. ()2.24|) and p.25p 
have the form Or = '4'i^'4^2 and the above Ward identities become 

dfj.Vfj,{x) ^ {mi - 1712) S{x) , (3.3) 

dM^) = (mi+m2) P{x) . (3.4) 

Beyond classical level, these Ward identities are to be understood as statements 
between operators. More precisely, they are insertions in expectation values of multi- 
local operators 0{xi, . . . ,Xn), consisting of a product of quark and gluon fields at 
different space-time points (xi ^ a;2 7^ • • • 7^ a;„). In the bare (lattice) theory these 
expectation values are defined in the path integral formalism as 

(0(xi,...,x-„)) = ^ J i'DU] [2?V^] [Vi'] 0{xi,...,xn) eM-SG-Sp] • (3.5) 

Lattice Ward identities are obtained in a way analogous to the continuum case. We 
recall that the basic trick is to consider, instead of the usual global SU{Np)l <E) 
SU{Nf)i^ chiral transformations, infinitesimal local vector and axial transformations 
of the fermionic fields; these are 



and 



Sip{x) 



S'ipix) 



A" 



ipix) 



S'ip{x) ~ -~i'ip{x) 



S'ipix) = iip{x) 



A" 



A° 

'5«l(a;)Y75 



(3.6) 



(3.7) 



respectively. The operator expectation value of eq. p.Sp is invariant under any change 
of the fermionic fields ^ and "05 since they are integration variables. This implies that 

^(0(.i,. ..,.„)) = 0, (3.8) 

where da"' may be either Say or Sa\, depending on the transformation under consid- 
eration (vector or axial). The above expression leads to the relation 



SO{xi,. ..,Xn) 
Sa'^{x) 



0{xi,. . .,Xn) 



SSf 



6a"-{x) 



(3.9) 



This is the quantum field theoretic expression leading to a Ward identity: for vector 
transformations Say, the variation of the fermionic action on the r.h.s. corresponds to 
the PCVC relation of eq. p.ip , while for axial transformations Say , it gives rise to the 
PCAC relation of eq. (|3.2p . The exact meaning of this statement and the subtleties 
that accompany it in the case of Wilson fcrmions is the subject of the rest of this 
section. 



^Note that unless otherwise stated, flavour indices run over / = !,■■■ , Np — 1. The singlet case 
(/ = 0), related to anomalous Ward identities, will not be discussed in these lectures. 
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3.1 Lattice Vector Ward identities 



It is easy to check that the global vector transformations corresponding to those of 
eq. ()3.6|) are a symmetry of the action ()2.2[) when quarks are degenerate in mass (i.e. 
Mo is proportional to the unit matrix). In common lore, "vector symmetry is a lattice 
symmetry". From the local transformations of eq. p. 61) the following vector Ward 



identity can be derived ( [Karsten and Smit, 1981[ ): 



./60{xi, 



0) 



^{x) 0(xi,...,.T„)) , (3.10) 



where aV^f{x) — {f{x) — f{x — fi)) is an asymmetric lattice derivative, and 



A" 

V'(2;)(7m ~ '^)U^i{x)—'il^{x + afi) + 
^{x + aA)(7A. + 



(3.11) 



is a point-split vector current. With x ^ xi, . . . ,x„ and in the limit of degenerate 
bare quark masses (Mq oc /), eq. p.lOp leads to the conservation of the point-split 
lattice vector current, V^V"^(x) = 0. The conservation of the standard local vector 

current V^(a;) = 'ip{x)—j^ip{x) on the lattice is somewhat more intricate, as we will 
see below. 

Without loss of generality, we again simplify matters by combining two versions 
of eq. (j3.10p . one with flavour index a ~ 1 and one with a = 2, into a single Ward 
identity, involving the raising Gell-Mann matrix A"*" of SU{Np)v and the two flavours 
ipi and ip2 (with corresponding bare masses moi and mo2)- The vector current and 
scalar density are then bilinear operators of the form Or{x) = ilJiTip2- We also chose 
the specific operator 0{xi, X2) = tpiixi)^p2ix2)- In terms of these quantities the lattice 
vector Ward identity (|3.10p becomes 



V^(^i'i{xi) V^{x) -02(a;2)y = [ rnai - mo2 ] (^ipi{xi) S{x) 'ip2{x2) 
- S{x2 - x) {ipi{xi)-ipi{x2)) - 6{xi - x) {ip2{xi)-ip2{x2)) ■ (3.12) 
Recalling the definition (|2.27p of the vertex functions Go, this is written as 

^^xG~ {xi -~x,X2- a;;moi,mo2) = (toqi - mo2) Gs {xi -x,X2- x;moi,mo2) 



+ d{x2 - x)Si (xi - X2; moi) - S{xi - x) S2 {xi - X2; TO02) 



(3.13) 



This expression relates bare quantities which are divergent in the continuum limit. Let 
us first work in the degenerate mass limit (moi = '7102), in which the first term on the 
r.h.s. vanishes. A glance at eq. ()2.18|) shows that multiplying the resulting expression 
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by Z^, renders the r.h.s. finite. Thus the l.h.s., which is now Z^^^^V'^G'i^, is also 
finite and does not require any further renormahzation. Consequently, the point-split 
conserved vector current satisfies the vector Ward identity with its normalization 

given by the trivial factor 

Zy = l . (3.14) 



This is the no-renormalization theorem on the lattice with Wilson fermions ( Karsten and Smit, 1981 
|Bochicchio et al., 1985[ ). It is an exact analogue of the familiar one in continuum QCD. 
The only difference is that the lattice vector current satisfying it is not the familiar 
one = "017/1 '02, but the point-split one defined in eq. p.ll[) . This result guarantees 
a proper definition of the vector charge and the validity of current algebra. 

We will now show that the correlation functions of the local vector current (x) = 
01 (2^)7^1 "02 (2;) differ from those of the conserved current by finite contributions, which 
vanish in the continuum limit. The argument goes as follows: first we express the 
conserved current as 



0i(a;)(7M - 1) M^) + Mx) (7m + ^)Mx) 
V^{x) + aA^{x) , (3.15) 



where we have used the lattice asymmetric covariant derivatives of eqs. (j2.3|) and 
(|2.5p . The second term on the r.h.s. of eq. p.lSp is a dimension-4 operator A^. For 
dcfinitiveness, we now consider the vertex function Gy{p); from eq. p.lSp it follows 
that 

G^ip) = Gvip) + aG^ip) , (3.16) 

which we have shown to be finite, up to an overall quark field renormahzation Z^; cf. 
eqs. (|2.34p and (|3.14p . The term aG^ip) vanishes at tree-level in the continuum limit. 
Beyond tree-level however, this term contributes, due to the power divergence induced 
by mixing with a lower dimensional operator. To see this, write the renormalized 
dimension-4 operator as 

r . „ r . c(gn, am) 1 

[A^Jr = Za [a^ + — . (3.17) 

This renormahzation pattern is dictated by dimensional arguments and the quantum 
numbers of A^. The dimensionless renormahzation constant Z\ is at worst logarithmi- 
cally divergent. The dimensionless mixing coefficient c{go,am) depends on the gauge 
coupling and the bare quark masses (collectively denoted by am). It has been shown 
that it cannot depend on a renormahzation scale a/i ( [Testa, 1998p . From the point of 
view of dimensional analysis, c{go,am) could contain logarithms of the form In(am). 
However, the absence of such a logarithmic dependence has been explicitly shown for 
the axial current at all orders in perturbation theory; the situation is analogous for 
the vector current ( [Curci, 1986[ ). This is not unexpected, as ln(am)-terms would lead 
to an ill-defined chiral limit at fixed lattice spacing. Thus we are left with a regular 
dependence of c{go, am) on the quark mass (i.e. positive powers of am), which may be 
dropped in a mass independent scheme leaving us with c{go). 
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From these considerations we deduce that, up to regular 0{a) terms (i.e. discretiza- 
tion effects), eq. p.l6p may be written as 

Gy{p) = Gv{p) + [Ga(p)]r - c(.go) Gv{p) 

= [1 - c{go)]Gv{p) + --- . (3.18) 

The term proportional to [Ga]r has been dropped in tlie last expression, being pro- 
portional to the lattice spacing (the renormalized correlation is finite by construction, 
while Za is at most logarithmically divergent). Consequently, the local current has 
a finite normalization: 

Zvigo) ^ l-c(.go) ^1 ■ (3.19) 

Note that eq. (|3.13[) (and any other vector Ward identity) can be expressed in terms 
of the local current, by substituting with [Zy V^]. 

Next we go back to eq. (|3.13p and, keeping the two masses distinct, integrate it 
w.r.t. a;; the surface trerm on the l.h.s. vanishes and we obtain 

(mo2 - moi) J d'^^ {xi - x,X2- a;;moi,TOo2) = 

Si{xi - X2;moi) - S2{xi ~ X2;mo2) ■ (3.20) 

Again, the r.h.s. of above expression becomes finite (renormalized) once we multiply 
through by Z^. This means that, up to this quark field renormalization, the l.h.s. 
is also finite. In other words, the product of the bare quark mass difference and the 
integrated scalar operator / d'^xS{x) is a scale independent (or renormalization group 
invariant) quantity. As there is no operator of dimension c? < 3 that could mix with 
S and vanish identically under the integral J d'^x, also the product of the bare quark 
mass difference and the (unintegrated) scalar density S is scale independent. With the 
difference of renormalized quark masses formally given in eq. (j2.16p and the operator 
renormalization given by eq. (j2.32p . this means that 

Zsigo,afi) ^ Z-^{go,afi) . (3.21) 

Thus, vector Ward identities imply an exact relation between the quark mass renor- 
malization Zjn and that of the scalar density. In other words, once a renormalization 
condition is imposed on the quark mass in order to determine Zm, the scalar density 
renormalization constant Zs is also known (or vice versa) . A corollary is that the mass 
anomalous dimension is the opposite of that of the scalar operator (7,„ = —7s). 

3.2 Lattice Axial Ward identities 

Far less obvious are the consequences of axial transformations with Wilson fermions. 
This is because even at vanishing quark masses (i.e. Mq = 0), the Wilson term of 
the lattice action p.2|) is not invariant under the global version of the axial trans- 
formations p.7|) . However, by imposing suitable renormalization conditions, PC AC 
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is recovered in the continuum limit ( [Karsten and Smit, 1981[|Bochicchio et ai, 1985[ 
[Testa, 1998[ ). The axial WI, namely eq. p.9p for the axial transformations p.7p . is 



■ ; ) 



-a'^(^(x-)|y,Afo|75^(a;) 0{xi, . . . ,Xn)) - a''(x^{x) 0{xi, . . . , x^)) , (3.22) 
where ^Ji(x) is a bilinear point-split axial current given by 

^"^i^) - I [^ix)wU^{x)Y^{x + afi) + 4,{x + a^)7^75t/t(x)y7/.(x-)] . (3.23) 

The term in the above Ward identity is the variation of the Wilson term under 
axial transformations: 



^ ' 2 



i^{x)—-izV ^{x)+4j{x)V Y75V'(a;) 



(3.24) 



where 



V i'ix) = ^ [Ufj,{x)^p{x + apL) + C/^(x - ajl)%}j{x - a/i) - 2il>{x)\ 
^2 



Since cannot be cast in the form of a four- divergence, it cannot be absorbed in 
a redefinition of the axial current A^. This is of course a corollary of the fact that 
the Wilson term in the fermionic action breaks chiral symmetry. It is a dimension- 
4 operator which, in the naive continuum limit vanishes, being of the form X" = 
aOg, with Og a dimension-5 operator. However, when inserted in correlation functions 
such as (|3.22p . Of will generate power divergences which will cancel the prefactor 
a. It is therefore mandatory to understand its renormalization properties. Its mixing 
with other dimension-5 operators may be ignored, since the corresponding mixing 
coefficients are at most logarithmically divergent and their contribution to = aOf 
vanishes in the continuum limit. Mixing with lower dimensional operators is, on the 
other hand, relevant, as the corresponding coefficients are inverse powers of the lattice 
spacing ( [Testa, 1998^ . The only possibility for such operator mixing, respecting the 
symmetry properties of Og is: 



0^{x) ]r = Zs 



(3.26) 



The functional dependence of the mixing coefficients and AI is established by 
arguments, very similar to the ones of sect. 13.11 concerning Zy- These coefficients 
involve subtractions of operators of lower dimension than that of 0§, which im- 
plies that they are independent of any renormalization scale /x ( [Testa, 1998[ ). Any 



18 Lattice Ward identities 



dependence on aAI cannot be divergent (e.g. In(aAf)) because this would compro- 
mise the chiral hmit of [Of (a;)]H-insertions in renormahzed correlation functions. The 
absence of ln(aM)-dependence of has been explicitly proved in perturbation the- 
ory ( [Curci, 1986[ ). Moreover, any regular dependence of Z^^ on aM will vanish in a 
mass independent renormalization scheme and may therefore be considered a lattice 
artefact. We conclude that the dimensionless axial current coefficient only depends 
on the gauge coupling, Z^{go). Similarly, the dimension-1 mass coefficient has the 
functional form M = w{gQ, aM)/a, with w{go,aM) a regular function of aM. 
Using cq. p.26|) . wc substitute = a x Og in eq. (|3.22l) . obtaining 

- a4(^(x)|y,[Mo-I/]|75V'(a:) 0(xi, . . . , x„)) 
-«'(^|^0(a;i,...,x„)) . (3.27) 

In order to simplify the presentation, we again specify the operator 0{xi,X2) ~ 
ipi{xi)ip2{x2), and the Gell-Mann matrix of SU{Np)a, in place of A". The above 
expression becomes 

V^l^ipiixi) Z^A^,{x) ■02(a;2)^ = [ rrioi +^02 - mi - m2 ] ^-01(2^1) Pix) ■4'2{x2)^ 

- 5{x - X2) {■lkl{xi)'l[}l{x2)lb) - &{x - Xl) {■~1<i'4>2{xi)-li)2{x2)) (3.28) 
+ ^ (^1(2:1) [05(2:)]r Mx2)) , 

with mi , m2 the elements of the diagonal! matrix M. Note that these quantities depend 
on the bare coupling and all bare quark masses; i.e. rnf{gQ, mgi, ■ ■ ■ , mQNp). The last 
term on the r.h.s. is genuinely 0{a). Since it vanishes in the continuum limit, it may 
be safely dropped as a discretization effect. This leaves us with the axial Ward identity 

-x,X2- x; moi, mo2) (3.29) 

(mf + m^^^^) Gp(xi-2:,X2-x;moi,mo2) 
- (5(x2 - x) 5i (xi - X2;moi)75 - (5(xi - x) 75 ^2 (xi - 2:2; mo2) + 0{a) , 

where the PCAC (bare) quark mass is defined by 

y^PCAC ^ ^^^^ _ fnf{gQ,moi,- ■ ■ ,moNp) (/ = 1, 2, 3, • • • , iV^p) . (3.30) 

We may now define the chiral limit on the basis of Ward identity p.29p : first we take 
for simplicity al masses mo to be degenerate, which allows us to simply the notation 

^If the matrix A/ were not diagonal, quantum numbers like strangeness would be violated by the 
theory. 
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by writing my as rn(go, mo). Then the chiral hmit is defined as the value mcr of mo for 
which m((7o,"icr) — ^ci (i-e. m^*^^*-^ vanishes). A relation between the PCAC quark 
mass defined in eq. p.30|) and the subtracted bare mass m = niQ — is readily 
obtained by expanding the former around the chiral point mcr: 



PCAC _ _ ^ dm{m) 



m — mo — m(mc 



dm 



(mo - mcr) + 



= (mo -mcr) h - + •••1 . (3.31) 

In the chiral limit, the first term on the r.h.s. of eq. p.29p vanishes. Multiplying 
the resulting expression by renders the r.h.s. finite. Thus, also the l.h.s. is finite 
and does not require any further renormalization. Now all terms of eq. (|3.29p differ 
from their continuum expressions by C(a) discretization effects. The point-split axial 
current p.23|) with the normalization 

Z^^l , (3.32) 

satisfies the continuum axial Ward identity up to C(a). In practice it turns out to be 
more convenient to work with the lattice local axial current A^(x). We can show, in a 
fashion analogous to the case of the vector current (c.f. the power counting argument 
based on eqs. ()3. 15^ - 1)3. 19|) ). that A^(x) has a finite normalization constant Za- Thus 
we have [A°]r = lima^o[-^^^J!l] = linia^o[-^A^y ■ From now on, the combination 
Z^A'^^ will always be substituted by ZaA'^^ wherever it appears in a WI. Also analogous 
to the vector current case is the lack of mass dependence of these constants. We 
therefore have: 

ZM) ' ZA{gl) + 1 . (3.33) 

Note that the tree level value of both Z-^ and Za is unity. In perturbation theory 
Za 1 + ^i5o + ■ ■ ■ (and similarly for Zj^. Thus in the continuum limit (a — > 0, 
g^(a) — > 0) we have Za-,Z-^^ 1. 

From here on things proceed very much like the case of vector Ward identities. Inte- 
grating eq. (j3.29p over x eliminates the surface term with the axial current divergence, 
leaving us with 

(mf^^^ +771^'^'^^) y" d^x Gp {xi - x,X2 - x;moumo2) = 

5 (xi - 2:2; moi) 75 + (.Ti - X2; mo2) . (3.34) 

Upon multiplying through with the quark field renormalization constant Z.^, the r.h.s. 
becomes finite. Thus also the l.h.s. is finite, which means that the product of the 
PCAC quark mass times the pscudoscalar density is a renormalization group invariant 
quantity (i.e. it is scale independent). The anomalous dimension of the the pscudoscalar 
density is the opposite to that of the quark mass (7,„ — — 7pjfl. Given that the 

*This is not unexpected: anomalous dimensions are continuum quantities and since in the con- 
tinuum scalar and pseudoscalar operators belong to the same chiral multiplet, they have the same 
anomalous dimension 73 = ■yp. We have already derived from vector Ward identities that 7s = —-fm- 
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pseudoscalar operator P is multiplicatively renormalized by Zp, this suggests the 
following renormalization pattern for the PCAC quark mass: 

[ "i/(gR,Ai) ]r = Zp^{ga,an) m^^'^'^{go,moi, ■ ■ ■ moNj,) ■ (3.35) 

This expression amounts to a definition of the renormalized quark mass, in a scheme 
which fixes the renormalization of the pseudoscalar density Zp. Identifying this defi- 
nition of [to/]r with the generic expression (|2.13p gives a closed expression for rrif. 

In conclusion, lattice axial Ward identities lead to: (i) the definition of the chiral 
limit as the solution of the equation m^'^^^(mcr) = Tier; (ii) the normalization of 
the axial current by Za- The last requirement may raise the following question: 
the properly normalized matrix element, from which the pion decay constant f^^ is 
computed, is given by 

Za (ttIA^IO) = (MiAMO) + 0{a) . (3.36) 

Since we know that in the continuum limit Za 1, why bother to normalize the 
axial current by Za'? The answer provided by axial Ward identities like eq. p.29p (and 
eq. (|3.38p below) is that close to the continuum limit the normalized matrix element in 
the above expression differs from its continuum limit by discretization effects which are 
positive powers of the lattice spacing. Consequently, omitting Za (a regular function of 
the bare coupling cjq), implies that the bare matrix element converges to its continuum 
limit like gQ l/ln(a). This is much slower than the convergence of ZaA/^^, shown in 

eq. ([Xse]) . 

3.3 Hadronic Ward identities for Zy and Za 

So far we have seen how Ward identities fix the normalization of lattice partially 
conserved currents and the ratio of the scalar and pseudoscalar renormalization con- 
stants. We have mostly used vector and axial Ward identities based on the variation 
of the operator 0{xi,X2) = ip{xi)'ip{x2)^ commonly referred to as "Ward identities on 
quark states" . This terminology refers to the fact that at large time separations, say 

< x° < a;^, expressions like (|3TT3)) and (|3^ eventually produce PCVC and PCAC 
relations between matrix elements of quark states. These Ward identities can form a 
basis for the determination of the current normalizations Zy-, Za and the ratio Zs/Zp^ 
but this is usually not very practical (for instance they require gauge fixing, which is 
affected by the Gribov ambiguity etc.). It is therefore preferable to write down Ward 
identities involving correlation functions of gauge invariant operators, based on suit- 
able choices of 0{xi,X2)- These are called hadronic Ward identities because they give 
rise, in terms of the LSZ procedure, to relations between matrix elements of hadronic 
states. Here we will briefly describe some examples of Ward identities leading to non- 
perturbative computations of the scale independent parameters Zv-,Za- A similar 
discussion concerning ratios of scalar and pseudoscalar renormalization constants is 
the subject of the next section. 

The computation of Zy is based on the observation that, up to discretization 
effects, = ZyVfj, + 0{a). It is then straightforward to obtain, at fixed UV cutoff 
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(i.e. fixed go) estimates of Zy by comparing tlie insertion of the point-split vector 
current in a correlation function, to that of the local current. For example, equations 



^ /d^x ( V^^ix) V^{z) ) = Zvj^ fd'^i V^{x) Vt{z) ) , (3.37) 

k=l •' k=l •' 

d^x / d^y { P^{x) V^{z) P^y) ) = Zy f d^x / d'y { P%x) Vo\z) P^y) ) , 



are valid up to 0{a). Thus they may be solved for Zy, providing two independent 
estimates for this quantity, which differ by discretization effects. Strictly speaking, 
integration over space dimensions is not an essential feature, but it is useful in practice, 
as it averages out statistical fluctuations. In order to avoid singularities from contact 
terms, time-slices are kept distinct (x-° 7^ 7^ 2"). Note that the flavour superscripts 
a, 6, c must be chosen so as to give, through Wick contractions of the fcrmion fields, 
connected diagrams of valence quarks with non- vanishing expectation values. 

A non-perturbativc determination of Za is based on two axial Ward identities. We 
write down cq. p.27|) for the operator 0{y) = P^{y); with y ^ x the l.h.s. vanishes. 
As we have argued above, the term containing may also be dropped, being a 
discretization effect. We also work with the local, rather than the point-split current. 
Thus, up to 0{a), we have 



/J 



^j{x) 



[A/o - M] 
2 ' Za 



P\y) 



(3.38) 



Like before, flavour indices a, b must be appropriately chosen. The above Ward identity 
may be solved for the ratio [Mq — Mo\/Za; this is one of the most standard methods 
for the determination of the PCAC quark mass, up to the axial current normaliza- 
tion factor Za- At fixed gauge coupling and degenerate quark masses, the quantity 
[mo — mo\/ZA is computed at several values of the bare quark mass mo and subse- 
quently extrapolated to zero, in order to obtain a non-perturbative estimate of mcr; 
cf. eq. p.30p and the comment that follows it. Note that in order to increase the signal 
stability, the above Ward identity is usually integrated over all space coordinates of y 
(with x° ^ ?/°). 

The next step is to work out the axial Ward identity (|3.22p . with the operator 
choice 0[x,y) = Al{x)V;;[y) (for Np > 2): 



Za ( Aliz) Atix) F;(y) ) 



(3.39) 



= ( [i}{z) {iA^ {Mo A/)} 75 ^{z)] A':^{x) V;{y) ) + ( X'^(z) A^x) V^^{y) ) 
+ z r'"' S{z - x) { V,\x) y;(y) ) + 2 r'' 5{z - y) { Al{x) A^(y) ) , 



where again use of eq. (|3.26p has been made and the shorthand notation X = 
a[Of\-R,/Zz has been introduced. Once more, flavour superscripts a, &, c are chosen 
so as to give, upon Wick contraction, non-vanishing expectation values of connected 
diagrams of valence quarks. The correlation function with the X operator gives rise 
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to contact terms. Symmetry arguments impose that these contact terms must have 
the form 



( Alix) 



^H9o) r" 5{z~v){Al{x) Aliv)) 



(3.40) 



where the elhpsis stands for focahzed Schwinger terms, which wih vanish after inte- 
gration over z, to be performed befow (wc always keep x ^ y). For x^y ^ z the l.h.s. 
vanishes in the continuum hmit as discussed in section [3.21 

We now proceed as follows: (I) rewrite eq. (|3.39p . using eq. (|3.40p and expressing 
all bare quantities in terms of the renormalizcd ones; (II) recall that the product of 
the PCAC quark mass and the pseudoscalar density in the first term of the r.h.s. 
of cq. p.39p is renormalization group invariant (c.f. eq. p.35p ): (III) require that 
the renormalizcd quantities obey the corresponding continuum (nominal) axial Ward 
identity. Thus we obtain 



fci(5o) = l- 



Za 



Za 
Zv 



1 



(3.41) 



and the Ward identity in terms of properly normalized currents is: 

A° M-Mn 



Za Zv 



{ [w^^A;{z)-i}iz){-,^^}j,i;{z)] Alix) F;(j/) ) = 



+ zr'"'Siz-x)Zl{V,'{x)V;iy)) 
+ 1 rU{z-y) Zl{Al{x) A^(y)) 



(3.42) 



Performing the integration over z kills off the first term on the l.h.slf|. We also integrate 
over X in order to improve the signal to noise ratio in practical computations (recall 
that X ^ y \& necessary in order to eliminate Schwinger terms; thus x^ ^ y'^): 



d^z I d^x ( [^Piz) {iA^^^^^^}75^(^)] Alix) F;(2/)) 



(3.43) 



' |r J { K'(^) V^iy) ) - ^ ^ I d^x ( AUx) ^^(y) ) 

Since Zv is known from eqs. p.37p and [M — Mo]/Za is known from eq. p.38p . the 
above Ward identity may be solved for Za- For Np = 2, this determination of Za is 
clearly not viable. A method for obtaining Za for any Np > 2 is based on an axial 
Ward identity with quark external states ( Martinelli et al., 19931. 



^ Note that eq. I|3.43| l is only valid away from the chiral limit, where the integral over z of the 
total divergence W^^lii^) vanishes. At zero quark mass, the term containing the total divergence of 
the axial current contributes, because of the presence of massless Goldstone bosons. 
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3.4 Hadronic Ward identity for the ratio Zs/Zp 

A further application of hadronic lattice Ward identities concerns the renormaliza- 
tion of the scalar and pseudoscalar densities. In a regularization which respects chiral 
symmetry, these operators, belonging to the same chiral multiplet, are renormalizcd 
by the same parameters Zs = Zp. This is not the case of Wilson fermions, which 
break chiral symmetry. Following a line of reasoning similar to the one of the previous 
section, applied to the operator 0{x,y) ~ S^{x)P^{y) (with g ^ h), wc obtain the 
Ward identity (for Np > 2) 



df^^ I d^x ( P\x) P\y) ) + d/« / d^^ { S^ix) S\y) ) 



ZaZs J Zj^Zp 

where the d's are defined in eq. ()2.7|) . Once Za has been determined through, say, 
eq. p.43p . the above Ward identity may be used to compute the ratio Zs/Zp. 

Note that Ward identities such as p.43p and p.44p may be regarded as explicit 
demonstrations of the fact that Zy,Za and Zs/Zp are scale independent functions 
of the bare gauge coupling: upon solving them for these Z-factors, we obtain solu- 
tions which are combinations of bare correlation functions. These bear no dependence 
on renormalization schemes and scales and remain finite in the continuum limit. At 
fixed gauge coupling (i.e. fixed UV cutoff), they provide non-perturbative estimates of 
Zv,Za and Zs/Zp. 

Practical simulations of these Ward identities are often performed on lattices with 
periodic and/or antiperiodic boundary conditions for the gauge and fermion fields. 
For a given gauge coupling, simulations are carried out at several non-zero quark 
masses and the resulting Z-factors are extrapolated to the chiral limit. Alternatively, 
Schrodinger functional (Dirichlet) boundary conditions and suitably chosen correlation 
functions enable us to perform computations of Zv,Za and Zs/Zp directly in the 
chiral limil[^. Possible differences among numerical results for the same Z-factors, 
obtained from different formulations and/or Ward identities, are due to discretization 
effects and provide an estimate of this source of systematic error. 

A determination of the ratio Zs/Zp, valid also for Np ~ 2, is based on the expres- 
sion 

[ ]r, - [ TO2 ]r = Zg^[moi - mo2] 

= Zpi[mPCAC _ ^PCAC] (3 45) 

The first line is obtained from eqs. (|2.16p and p.2ip . while the second one from 
eq. p.35|) . The ratio of bare to PCAC mass differences gives an estimate of Zs/Zp. 
Such a computation must be carried out for pairs of non-degenerate quark masses. 

''For scale dependent renormalization parameters, such as Zp, renormalization group running 
is naturally performed non-perturbatively in the Schrodinger functional scheme, for a large range 
of renormalization scales. As this topic is beyond the scope of the present lectures, the reader is 
advised to consult the literature on non-pertrubative renormalization and the Schrodinger Func- 
tional ULuscher, 1998l|Sommer, 1997||Weisz, 2010[ l. 
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Subsequently results are to be extrapolated to the chiral limit. To the best of our 
knowledge, this method has never been applied in practical computations. 

3.5 Singlet scalar and pseudoscalar operators 

It is sometimes convenient to work with Ward identities involving matrix elements 
of gauge invariant operators between hadronic states, obtained from Ward identities 
between correlations functions such as eq. p.39|) . and the use of the standard LSZ 
procedure. The reader should have no difficulty performing these steps. We focus on 
the following Ward identity, based on the axial variation of the pseudoscalar density 
amongst hadronic states \hi) \h2): 



(h^inx'ix) p^ 



+ i 



hi 



2 [too - 
dPs{0) 



to] {hi\T[P^ix) P3{ 



(3.46) 



The flavour index / = 1, • • • , A^J. — 1 characterizes non-singlet currents, pseudoscalar 
densities etc., while g = 0, • • • ,N'p — 1 also covers the case of the flavour singlet 
pseudoscalar density P*^. For simplicity we assume exact flavour symmetry; i.e. A/o = 
toq/ and thus {A-'^, Mq] ~ 2mQ\f . The axial variation of the pseudoscalar density is 



5PS{Q) 



(3.47) 



The axial current normalization Za and the mass subtraction to in eq. p.46p are 
generated in standard fashion by the counter-terms of the operator O5 = /a; cf. 
eq. (|3.26p . This is not, however the end of the story as we now allow for the case in which 
the space-time point x comes close to the origin (a; ~ 0). In this case, the insertion of 

X in off shell green functions does not vanish as we approach the continuum limit. 
Rather, it generates local contact terms which are constrained, by flavour symmetry, 
to have the form: 

{hi\T[X^ix) P^m \h2) = - Six) \ cligl) d^'" (hilS'^mh^) (3.48) 



C3(g^)'5^y^ {hiisfmh. 



Putting it all together, and going over to the chiral limit, yields the Ward identity 
{hi\T[ Za Alix) P^O)] \h,) = 5{x) { [l + ci(5o')] d^'" {hi\SHO) Ih^) 

-C2igl)] {hi\S"iO)\h2) (3.49) 

cMg'o)] {hi\SHO) \h2) } . 

For a non-singlet pseudoscalar density P^ (i.e. for g 0) and for f ^ g, the 
last terms in the r.h.s. of the above expression vanish. Since P^ is multiplicatively 
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rcnormalizable, it is adequate to multiply the above expression by Zp to render it 
finite. This results to the following relation between the renormalization factor of the 
non-singlet scalar density and Zp: 



Zs = [l + c,{gl)] Zp 



(3.50) 



If, on the the hand, f = g ^0, the second term on the r.h.s. of eq. p.49p survives and 
we obtain the relation: 



Zgo = [l + C2{gl)] Zf 



(3.51) 



Finally we consider the case in which is a singlet pseudoscalar density (i.e. g = 0), 
for which only the last term on the r.h.s. of eq. (|3.49p survives. Multiplying through 
by the factor Zpo wc obtain the relation 



Zs 



[1 + C3{gl)] Zpo 



(3.52) 



Thus it has been established ( [Bochicchio et ai, 1985|Rossi et ai, 2003[ ) that the [Np, Npji 
[NpjNp] representation of the chiral group SU{Np)l (X) SU{Np)p is formed by the 
rescaled operators 



pf 



1 + Ci 



p.. 



[1 + ci] 

[1 + C2] , 



(3.53) 
(3.54) 



which renormalize by a common factor Zp. Conversely, the renormalizcd operators 
pI = ZpPf, = ZgPf (for / 7^ 0) and P^ = ZpoPf have a common anomalous 
dimension (arising from, say, Zp) but renormalization constants Zp, Zg and Zpa which 
differ by finite normalization factors. These factors tend to unity in the continuum 
limit. Their presence is yet another consequence of chiral symmetry breaking by the 
Wilson term. Analogous statements are true for the renormalization factor Zga of the 
singlet scalar density , which however is also subject to power subtractions, as we 
will see below. 

Another important result (Maiani et ai, 1987 Rossi et ai, 20031 is the relation 



between the ratio Zp/Zgo and the ratio of PCAC to subtracted quark masses. The first 
step towards its derivation consists in a Ward identity analogous to that of cq. p.46|) . 
but with S° in place of 

{h,\r[ZA M{y) S°{x)] \h2) - 2 K - m] (/ii|r[P^(y) S^ix)] \h2) 



= {hi\T[x\y) S°{x)] \h2) + I (h. 



(550(a;) 



= -5{x-y) [l + tigD] 



{hi\P^{v)\h2) 



(3.55) 



The last equation is the result of performing the functional derivative 5S'^{x)/5a'^j^(ij) 



and establishing that the contact term, arising when X (y) is in the vicinity of S'^{x) 
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is proportional to P-^ , with t(5o) the proportionality factor. Multiplying through by 
Zgo gives the renormalized Ward identity 



ZsaV; {h,\T[ZAAl{y) S"{x)]\h2) - 



2 [mo - m] Zso {hi\T[P^ {y) j/is) 

(3.56) 



with the identification Zp [1 + t{g^)]Zso. 

The next step consists in differentiating the axial Ward identity between hadronic 
states with respect to the bare quark mass: 



d 



■VC {h,\ZAAUy) \h2 



d 



dmo ^ ^ dmo 

The mass derivative on the r.h.s. gives 
d 



2 [mo - m] {hi\pf{y)\h2 



(3.57) 



dmo 



2 [mo - m] {h,\pf{y)\h2 



1 



dm 



dmo 
2 [mo — m 



{hi\pfiy)\h2 
d 



dmr 



{h,\pf{y)\h2) . (3.58) 



Now differentiating a Green's function with respect to the quark mass mo amounts to 
inserting the operator J d'^xip{x)ip{x) = \/2Nf J d'^xS'^{x). Performing this insertion 
on the l.h.s. of eq. p.57p and the last term on the r.h.s. of eq. p.58p results to the 
following identity 



{h,\ZAT Aiiy) / d^xS'^ix) \h2) 



dm 
dmo 



{hi\PHy)\h2 



2 [mp — m ]{hi\T\p^{y) d^xS°ix)] \h2) . (3.59) 



Multiplying both sides of the above by Zgo gives a finite expression. Upon comparing 

it with eq. p.56p . and taking the chiral limit mo — > mcr, we obtain a relation between 

rcnormalization factors 

Zp _ dm 

Zso 



= 1 



dmo 



(3.60) 



Compared with eqs. p.30p and p.3ip . this gives a relation between the PCAC and the 
subtracted quark masses: 

,^PCAC 



lim 



Zp_ 

Zso 



(3.61) 



-mo - mcr- 

Combining this expression with the renormalized quark mass definition p.35p . we 
obtain 



rriR 



Zaa [ mo 



(3.62) 



which, upon comparison with eqs. (|2.14p and (|2.15p leads to 

Zso = Zj (3.63) 

We see that the above expression, as well as eq. p.21|) . relate the singlet and non-siglet 
mass renormalization constants to the corresponding scalar density rcnormalization 
constants. 
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In perturbation theory, the difference between singlet and non-singlet operators 
arises from extra diagrams in the singlet case, involving fermion loops. As these drop 
out in the quenched approximation, the quenched renormalization constants obey 
Zg = Zga and Zp = Zpo. 

3.6 The chiral condensate 

Our last topic related to lattice Ward identities is the derivation of the proper definition 
of the chiral condensate with Wilson fermions ( [Bochicchio et ai, 1985[ ). The starting 
point is the Ward identity (|3.46p for the vacuum expectation values of the non singlet 
axial current Aj^ and the pseudoscalar density (i.e. /, .g 7^ and = I/12) = |0)): 

(0| T[Za Alix) P^O)] |0) - 2 [mo - m] (0| T[P^ (x) P^O)] |0) 
= (0| T[X^ix) PS(0)] |0) + S{x) (Oli'MO) • (3.64) 

Note that only the term proportional to S*" survives in the vacuum expectation value 
of the axial variation of the pseudoscalar density p.47p . Moreover the contact terms 
arising in the limit of vanishing lattice spacing from the insertion of X'^ (x) with P^(0) 
at x w are not exactly those of eq. (|3.48p : at first sight it appears that the only 
contact term surviving in eq. ()3.48|) . sandwiched between vacuum states is the one 
proportional to (0|S'''|0). But this is not the whole story, since we should also take 
into account terms which vanish in eq. ()3.48|) . when hadronic states \hi) 7^ I/12), but 
survive when these states are equal (e.g. the vacuum). The reader can easily convince 
himself that in the chiral limit the contact term structure is 

(0| T[X^{x) P9(0)] |0) = [-^ 6{x)^^^ + □5(.T)^iM)l ^ (3.65) 

iNp a'^ a J 

where the term proportional to (0| 5*^10) is subdominant (it may be considered as 
incorporated in bo). Away from the chiral limit we also have terms proportional to the 
mass (e.g. cx m^'^^^/a^) but these are subdominant compared to the cubic divergence 
and irrelevant for the present discussion. 

Combining the last two equations in the chiral limit we obtain 

(0| T[Za Alix) Pf(0)] |0) - nSix)^^^ 

= S{x) \{0\4>m + 1 . (3.66) 

Np I a-^ i 

The subtraction on the l.h.s. compensates a divergent term, arising in the expectation 
value (0| T[Al{x) P^{0)] |0) when x w 0. This term is proportional to d^S{x). Muhi- 
plying the Ward identity by Zp renders both sides finite. This leads to the following 
definition of the renormalized chiral condensate with Wilson fermions 



1 



(^^)r ^ Zp— (OIV^IO) + . (3.67) 



a3 
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Note that a chirally symmetric regularization would lead to the much simpler result 
(^'0)r = Z so {0\'ilJip\0) / N f . The more complicated renormalization pattern of the last 
expression is due to the loss of chiral symmetry by Wilson fcrmions. 

Useful information may be also be obtained from Ward identity p.64p . by inte- 
grating it over all space-time. We distinguish two cases. In the first case we take the 
chiral limit before performing the integration. Then the term proportional to (toq — rn) 
vanishes. The integrated V^'A^-term would also vanish, being the integral of a four- 
divergence, provided chiral symmetry were not broken in QCD (absence of Goldstone 
bosons). Since upon integration, the Schwinger 6i-term of eq. p.65p also vanishes, 
the integrated Ward identity p.66|) would then imply a vanishing chiral condensate. 
However, in QCD the symmetry is broken, and the presence of Goldstone bosons guar- 
antees a non-vanishing surface term upon integrating V^A^. Thus, also the r.h.s. of 
the Ward identity (i.e. the chiral condensate) is non-zero. 

The second case of interest consists in integrating Ward identity p.64p over all 
space-time, before going to the chiral limit. Now this integration will kill off the first 
term on the l.h.s. (it is an integral of the total derivative of the axial current in 
the presence of massive states) while the second term survives. Upon integration, 
the Schwinger 6i-term of eq. p.65|) also vanishes. Taking the chiral limit after the 
integration, we obtain 



Sf9 



N 



f 



{OH'^\Q) + = - lim 2[mo-m] I d^x{0\ r[P^ {x) P^{Q)] |0) . 

J mo-i-mcr J 

(3.68) 

From the above equation, we can derive two other expressions for the chiral con- 
densate, by inserting, in standard fashion, a complete set of states in the time-ordered 
product of pseudoscalar densities. The spatial integration / d^x projects zero- momentum 
states. Contributions from higher mass states vanish in the chiral limit, leaving us with 
a zero-momentum pion state |7r(0)). Upon performing the time integration J dx'^, we 
find 



1 

N~f 



mm + ^ 1 = hm l!?!L^ (0|F(0)|7r(0)) ' , (3.69) 



where is the mass of the pseudoscalar state |7r) and the flavour indices f,g have 
been suppressed for simplicity. 

Next recall that the definition of the pion decay constant is given by 

{0\[A^ix)]n\7:{p)) = i U Pf,cxp{ipx) . (3.70) 

We also know that 

ZaV^ {0\A^ix)\Trip)) = 2[mo-m](0|P(a;)|7r(p)) = - Uml cxpitpx) , (3.71) 

where the first equation is the axial Ward identity between states |0) and |7r), while 
the second one is obtained by derivation of eq. p.70|) . Taking the square of the second 
equation at p = and combining it with eq. p.69p we find 



1 

7V7 



(OIV^V^IO) + 1 = hm ^ . (3.72) 

o-^ J mo^iMcr 4(too — m) 
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Multiplying both sides by Zp yields the familiar Gell-Mann-Oakes-Renner relation 
dGell-Mann et al, 1968[ ): 

(^V)r = 1™,T^ • (3-73) 

mR->0 4mR 

Note that the non- vanishing of the chiral condensate in the last two equations implies 
the well known linear dependence of the pseudoscalar mass squared on the quark mass, 
close to the chiral limit. Provided that numerical simulations can be performed close 
to the chiral limit, eqs. (|3.69p an d p.72[) may be used for the computation of the chiral 
condensate with Wilson fermiontU 

This concludes our discussion of lattice Ward identities with Wilson fermions. A 
related subject is that of Ward identities of the singlet axial current, related to the 
question of C/(l)^ anomaly and the rj' mass. Although extremely important, this topic 
is beyond the scope of the present lectures. 



''Computing the chiral condensate directly form the trace of the quark propagator is not viable 
with Wilson fermions, due to the presence of the cubic divergence proportional to bo- 
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Momentum subtraction (MOM) 
schemes 



Wc have seen that Ward identities may be used in order to determine non-perturbatively 
the normahzation of partially conserved currents and finite ratios of rcnormalization 
constants of operators of the same chiral multiplet. These quantities would be iden- 
tically equal to unity if chiral symmetry were preserved by the lattice regularization. 
On the other hand, local operators like S"", P° and T"^, are subject to multiplicative 
rcnormalization by scale dependent rcnormalization parameters Zs, Zp and Zt, which 
must be fixed by a more or less arbitrary rcnormalization condition. We now present 
a rcnormalization scheme, known as the RI/MOM scheme, which is suitable for the 
non-perturbative evaluation of these parameters. 

Before presenting the RI/MOM scheme explicitly, we will discuss a few general 
features of momentum subtraction schemes. For dcfinitiveness, let us consider the 
dimension-3, multiplicatively renormalizable operator Op defined in eq. (j2.23p . Mim- 
icking what is usually done with operator rcnormalization in continuum momentum 
subtraction (MOM) schemes, we impose that a suitable renormalized vertex function 
between say, quark states, at a given momentum scale /x, be equal to its tree level 
value. For instance, 

"( p' I I P > 

(4.1) 

where \p), \p') are single quark states of four- momenta p,p' respecively. This is a 
familiar example of a momentum subtraction scheme (MOM-scheme) . Some important 
properties of such schemes are: 

• The rcnormalization condition, imposed on quark states, is gauge dependent. 
Gauge fixing is required and one typically opts for the Landau gauge. Care should 
be taken that no problems arise from such a choice. For instance, dependence 
of the results on Gribov copies of a given configuration ensemble is obviously 
undesirable. 

• MOM is a mass independent, infinite volume rcnormalization scheme. This means 
that in principle the condition (j4.1[) is written at infinite lattice volumes and 
vanishing quark masses. In practice numerical simulations are performed at large 
but finite volumes and non-zero quark masses. Results are then extrapolated to 
the chiral limit; their independence on the volume must be carefully checked. 

• The above renormalization condition is imposed in the chiral limit. By working 
at non-vanishing exceptional momenta, any infrared problems associated with 



^ = Zoia^i) {p'\0^\p )barc = {P'\0^\p )t. 
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vanishing quark masses are avoided ptzykson and Zuber, 1980] ) . 
• Once the bare matrix element in the above MOM condition is regularized on 
the lattice, the renormalization constant Zoiaji). computed non-perturbatively, 
is subject to systematic errors due to discretization effects which are typically 
0{afj,), ©(oAqcd) and C(am), for non-improved Wilson fcrmions. Effects which 
are 0{am) are extrapolated away by going to the chiral limit, while those which 
are 0{a^) and ©(aAqco) are entangled with the non-perturbative definition of 
Zo. 

The reliability of non-pcrturbativc Zq estimates, computed in a MOM scheme, de- 
pend on the existence of a so-called "renormalization window" for the renormalization 
scale /i, defined through the inequalities 

Aqcd < M < 0{a-^) . (4.2) 

The upper bound guarantees that 0{afi) discretization effects are under control (and 
so are ©(oAqcd) ones). The lower bound is imposed for two reasons. The first reson is 
that in many cases, related to operator weak matrix elements, the physical amplitude of 
interest A is expressed in terms of an operator product expansion (OPE), schematically 
written as: 

A = (/iHeffN) = {f\0{^^)\^)n , (4.3) 

where \i) and |/) are physical initial and final states. The scale dependent Wilson 
coefficient Cw comprises all short distance effects of the physical process under con- 
sideration. It is known (typically to NLO) in perturbation theory. The renormalizcd 
weak matrix clement 

{f\0{^l)\^)K = \iuiZ[go{a),a^i] (/|0[,9o(a)] |*)barc , (4.4) 

a— >0 

is a long distance quantity, computed non pcrturbatively. Its scale dependence cancels 
that of Cw (to, say NLO), so as to have a scale independent physical amplitude A. 
Since Cw is a perturbative quantity, it should be calculated at a scale /x, well above 
-'^QCD- But fi is also the scale at which Zq must be computed non-perturbatively, 
which leads to the lower bound of the renormalization window. 

The second reason for requiring that Aqcd fi is related to the existence of 
unwanted non-pcrturbativc effects, which affect the determination of some Zq^s at 
scales close to the infrared. This so-called "Goldstonc pole contamination" will be 
treated extensively below. 

Moreover, the UV cutoff must be well above the energy scales of the problem, 
such as Aqcd, fJ- (and any quark masses of active flavours). Finally, all these scales 
must also be higher than the lattice extension L, which acts as an infrared cutoff. Thus 
we must tune the bare coupling .go (a) so that 

L^^ < Aqcd < /x < a"^ . (4.5) 

With prcscnt-day computer resources, this hierarchy of scales is not always easy to sat- 
isfy in practice. One could then turn to finite-size scaling methods, similar to the ones 
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based on the Sclirodinger functional ( jWeisz, 2010[ ) . Although these methods have been 

described in general terms for the RI/MOM scheme some time ago ( [Donini et al., 1999 1, 

they have only been systematically formulated and developed very recently ( Arthur and Boyle, 2010[ ), 

with extremely promising preliminary results. 



4.1 The RI/MOM scheme 

A specific renormalization scheme, which is particularly well suited for non-perturbative 
computations of renormalization constants is the RI/MOM scheme. As the name be- 
trays, it is a momentum-subtraction scheme, similar to the one described above (the 
reason for the acronym "RI" is described in Appendix [9l The difference lies in the 
fact that the renormalization condition is not imposed on operator matrix elements 
of single-quark states, but on the amputated projected correlation function defined in 
eq. (|2.30[) . Since the renormalization pattern of this correlation function is given by 
eq. p.36|) . it immediately follows that the renormalization condition 



['^oif^^gRi'mn = 0) 



lim 

a-i-O 



Z^^{aiJ.,go) Zo{an,go) To(j},go,m) 



m 



> 



1 



(4.6) 

fixes the combination Z^^ Zq- Note that unity on the r.h.s. stands for the tree level 
value of the correlation function Tq , which has been constructed so as to satisfy this 
property. 

There are several ways to disentangle the two renormalization constants from the 
product Z^^ Zq- a conceptually straightforward method is based on the conservation 

of the vector current V^. We assumc0 that the RI/MOM condition for the conserved 
current is compatible with the Ward identity result of ea. (|3.14p . Zy — 1. This means 
that 



[ r^(/x,5R,mR = 0) ]r 



lim 

a->0 



9 9 



1 . (4.7) 



The quark field renormalization Z^, cancels in the ratio To{tJ-)/^Y{fJ')^ "which can be 
solved for Zq- Alternatively, one often uses the product ZyVf^ in place of V^, with Zy 
known from some Ward identity. 

Another method is to first compute Z^, from the RI/MOM condition of the quark 
propagator, based on cqs. (|2.19p and (|2.21l) : 



rE(A^) 



= lim 



m 



= 1 



(4.8) 







Once Z^, is known, eq. (j4.6p gives us Zq- This method is avoided in practice, because 
the definition of Fx; involves derivatives, which look rather ugly with discretized mo- 
menta. A slightly different renormalization condition is preferred for the fcrmion field 
renormalization: 



'^This assumption will be proved in section [4. 3 1 
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— Tr 
12 



= lim (Z'l,)-^ — Tr 



m 



> 



= 1 



(4.9) 



Operator renormalization constants Zq, computed with Z'^, rather than Z^, are said to 
be in the RI'/MOM scheme. In general we expect that Z'^^ = Z.^, + ©((/q). However, in 
the Landau gauge it is known that Z'^ = Z,p + 0{gQ). As this is the gauge of preference, 
differences between RI/MOM and RI'/MOM results are expected to be small. 

For the quark mass, an RI/MOM renormalization condition, based on cqs. (|2.20p 
and would be: 



r,„(/i) 



= lim 



>0 L 



Zj,^{an,go) Z,„^(a^,5o) ^m.{p,9o,'m) 



9 ? 

m ^ 



= 1 



(4.10) 



In numerical simulations, this is also inconvenient, as Tm requires derivation of w.r.t 
the mass. The following condition is used in practice: 



rm(M) ]r 



lim 

a->0 



Z.J{an,go) Z,^^{a^i,go) T'^{p,go,m) 



P - 
m 



1 



(4.11) 



with 



1 



12m 



Tr[5-i(p;i 



(4.12) 



Note that the renormalization pattern of FJ^^ is analogous to that of eq. (|2.22p for F„i. 
The two conditions are equivalent, as differentiation of eq. (|4.1ip w.r.t. the quark mass 
m results to eq. (14.101 ) 

In practical simulations, the RI/MOM' conditions are solved for Z'^ and Zq at 
fixed bare coupling. The lattices are big with periodic and/or antiperiodic boundary 
conditions. At each go, one has to work at several small quark masses and then ex- 
trapolate the results to the chiral limit. The whole procedure is repeated for several 
lattice momenta ap in order to empirically establish a range of renormalization scales 
ap = a^, within the renormalization window of eq. (j4.5l) . 

4.2 Goldstone pole contamination 

It is clear from the previous section that RI/MOM is a mass independent renormal- 
ization scheme, i.e. the renormalization conditions are imposed in the chiral limit. 
For field theories with vertices with degree four (such as QCD), renormalized at 
some fixed Euclidean point, it is well known ptzykson and Zuber, 1980| , that if all 
masses of a Feynman diagram go to zero, infrared singularities may appear, un- 
less the external momenta are non-exceptiona0. This is not our case: the renor- 
malization condition (j4.6|) is imposed on the amputatcd-projected Green function 
Toip), derived from the momentum-space correlation function Go{pi,P2) of eq. (|2.28p 



correlation function has non-exceptional momenta if no partial sum of the incoming momenta 
Pi vanishes. 



34 Momentum subtraction (MOM) schemes 



for pi = P2- In other words, the operator O carries zero momentum (q = pi — 
P2 = 0) and quark-field external momenta are exceptional. In order to investigate 
which correlation functions are subject to dangerous infrared singularities, wc apply 
the LSZ reduction formula ( [Peskin and Schroeder, 1995[ ) to the correlation functions 
Go{pi,P2) dPapinutto, 200 Returning momentarily to Minkowski space-time, we 
obtain for the pseudoscalar operator P: 



Gp{pi,P2) 



d^X2 cM^P2X2) {0\T[iPiO)4,{x2)Mq)) 



^6l(g") 
q-^ — 



(7r(g)|[P(0)]R|0) 



(0|[P(0)]R|7r(-g)) 4^ {n{-q)\T[mi'ix2)m 



(4.13) 



where the ellipsis indicates terms without pion poles; isospin indices arc implicillf|. 

Combining cqs. p.7ip and p.73p wc obtain for the vacuum-to-pion pseudoscalar 
matrix element 



(0|[P(0)]R|7r(p)) 



■ cyupiipx) 



(4.14) 



Since (V'V')r ^ (spontaneous symmetry breaking), the above matrix clement survives 
the chiral limit. This last result, applied to cq. (j4.13|) in the limit q^ 0, implies that 
the correlation function Gp{p,p) (and thus also rp(p)) have a 1/mJ pole. The presence 
of this so-called Goldstonc pole means that the RI/MOM rcnormalization condition, 
which strictly speaking is valid in the chiral limit, can only be applied to Zp with 
some precaution (see below for details). 

The situation is less dramatic in the case of the axial current, for which the LSZ 
reduction formula is 



G^(Pl,P2) 



d^X2 exp{ip2X2) {0\T[i;{0)i}{x2)]\Tr{q)) 



ie{q") 
q^ — 777^ 



(7r(-q)|T[V^(0)^(a;2)]|0) 



(^(g-)|[A^(0)]R|0) 



(4.15) 



Since the relevant matrix element behaves like (0| [^^(0)]R|7r(q)) ~ q^ (cf. cq (|3.70p ) 
in the limit — > the Goldstonc pole contribution of G'^{p,p) vanishes. However, 
as wc will sec below, the pole of G'^{p,p) still plays an important role in axial Ward 
identities, used for the determination of Za- Other correlation functions such as Gy, 
Gs and Gt (as well as Fy, Tg and Tt) have analogous LSZ reduction formulae. In 
these cases however, the interpolating operators Vq, S and T^^i, have Lorcntz properties 
which ensure that the corresponding intermediate single-particle state is not a pion. 



^Contrary to what we have done so far, we find it convenient to identify xi with the origin (i.e. 
x\ = 0) in the correlation function Gp{x\ — x,X2 — x) of eq. I|2.27|l . Then the Fourier transform, 
analogous to that of eq. 112.271 1. is performed w.r.t to q (conjugate variable of x) and p2 (conjugate 
variable of X2), in order to obtain G(pi,p2), with p2 = pi — q. 
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As the mass of these particles does not vanish in the chiral hmit, these correlation 
functions are free of Goldstone pole singulariticfl 

It is important to realize that the Goldstone pole term of Gp{p,p) vanishes in the 
limit of large external momentum — >■ oo. This has been discussed in general terms in 
the original RI/MOM paper ( [Martinelli et ai, 1995[ ), where the following correlation 
function was considered: 

Foip) = J d^x cM-ipx) Oriy)) , (4.16) 

with Or{y) defined in cq. ()2.23|) . In the limit of large p^, the dominant contribu- 
tion to Fo{p) comes from regions of integration where the integrand is singular; 
i.e. a; ~ and a; ~ y. Based on the OPE and dimensional arguments, it is shown 



that (Martinelh et al, 1995) 



,M ^ + * A.(0, , (4.17) 

where cr, 7r, dr and Sr can be calculated in perturbation theory, while 

Ar(0) = / d^y (Or(0) Or(y)) (4.18) 



is a non-perturbative quantity. This means that, in the large limit, the perturbative 
contribution of Fo{p) (i.e. the first term, proportional to cr) dominates by one power 
of p^ over the non-perturbative one (i.e. the term proportional to dr). This is why 
Foip) is infrared safe in perturbation theory, when p^ is large. In other words, in the 
large external momentum limit, the contamination by the Goldstone pole residing in 
Ar(0), as well as other infrared chiral effects, vanishes. At finite momenta however, 
the Green function Fp, defined in eq. (|2.30p . is not regular in the chiral limit. By 
requiring that Aqcd << fJ- (the lower bound of the renormalization window (j4.5l) ) 
one is attempting to reduce, in practical simulations, the contamination due infrared 
chiral effects, such as the Goldstone pole. 

4.3 RI/MOM scheme and Ward Identities 

Once the wave function renormalization is computed, the RI/MOM renormaliza- 
tion condition (j4.6|) fixes the renormalization parameter Zq; examples are Zs, Zp^ Zt 
but also Zv and Za- However, Zy, Za and the ratio of Zs/Zp do not depend on any 
renormalization scheme (or renormalization scale), as they are fixed by Ward identities. 
It is then important that their determination through the RI/MOM renormalization 
condition be compatible with such identiticfU. In this section we will establish that 

^Intermediate two-pion states arc of course allowed. These however give rise not to poles, but to 
branch cuts with less severe, logarithmic singularities in the chiral limit. 

^Strictly speaking, it is of course conceivable that one renormalizes, say S and P, without taking 
Ward identities into account. This is acceptable from the point of view of renormalization (i.e. removal 
of divergences) but distorts, by terms finite in the bare coupling go, the recovery of symmetry as the 
regularization is removed. 
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the RI/MOM determination of Zs/Zp agrees with that from Ward identities only in 
the hmit of large renormalization scales fi » Aqcd, where Goldstone pole contami- 
nations die ofo. On the other hand, the RI/MOM determination of Zy is compatible 
with Ward identities at all momentum scales. Finally, Goldstone poles do not affect 
the RI/MOM determination of Za, provided the chiral limit is approached with care. 

Let us start from the vector Ward identity p.l3|) . with replaced by ZyVf^. We 
Fourier-transform and amputate eq. p.l3|) . obtaining (at small ag^) 



, 9 ^ 



S ^[p + -;moi 



(4.19) 



We next go through the following steps: (i) derive the above expression w.r.t. g^; (ii) 
let moi = mo2, so as to dispose of the term with the scalar vertex function A5; (iii) 
trace the resulting equation with 7^, (iv) take the limit g^ — )■ 0. This gives 



ZyTvip) + Z\, 



lim — Tr 



9k- 



^0 48 



7p- 



9Ay(p + g/2,p-g/2) 



= rs(p) 



(4.20) 



The LSZ reduction formula for the vector current is analogous to eq. ()4.15|) . but the 
denominator does not involve a Goldstone pole m^r. This implies that, upon taking its 
derivative w.r.t. gp, and in the limit — > 0, the second term on the l.h.s. of the above 
Ward identity vanishes (there are no infrared singularities). Thus we obtain 



Zy Vy{p) = rs(p) 



(4.21) 



Solved for Zy this is yet another Ward identity determination of the vector current 
normalization. On the other hand, if the above expression is multiplied by Z^^ ^ de- 
termined in the RI/MOM scheme through cq. (|4.8p . we obtain at momenta p ~ fi 



Z.J Zy Ty{n) 



1 



(4.22) 



which is the RI/MOM condition for Zy. We have shown that for the vector current 
the RI/MOM scheme is consistent with Ward identities. 

The situation is not quite the same for axial Ward identities. Starting from eq. (|3.29p . 
written for degenerate quark masses, we perform the same steps which led from 
eq. dHH) to eq. and obtain 



Za ^a{p) 



Za lim ^Tr 

q^^O 48 



2m 



757p- 



aAi^(p + g/2,p-q/2) 



PCAC 



48 



■ lim Tr 



757p- 



dqp 

dAp{p + q/2,p^ 



(4.23) 



9/2) 



dqp 



rs(p) 



®For simplicity we limit our discussion to the renormalization properties of quark bilinear op- 
erators. Analogous conclusions may be drawn for other more complicated cases. For example, the 
scale independent mixing coefficients of the fou r-fermion operators, such as the ones invo lved in 
K — ^ TTTT decays, are also fixed by Wa rd identities l|Bochicchio et al, 1985"]|Donini et al., 1999[ |. When 
computed in the RI/MOM scheme | |Donini et al., 1999| l, they are also subject to Goldstone pole 
contamination l |Becirevic et al., 2004[ l. 
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The LSZ reduction formulae (i.e. eq. (|4.15|) for the axial current and eq. (|4.13|) for the 
pseudoscalar density) may now be used in order to study the Goldstonc pole contri- 
butions to the second and third term on the l.h.s.. It can be easily shown that, away 
from the chiral limit, as — >■ 0, these terms vanish. Approaching subsequently the 
chiral limit, just like in the vector current case, the axial Ward identity determination 
of Za is equivalent to the RI/MOM condition. If however, the chiral limit m^'^^'^ — )- 
is taken before — 0, the third term vanishes but the second one survives. This 
term is essential to the saturation of the Ward identity, in the presence of a massless 
Goldstone boson (i.e. any divergences due to the presence of massless pions in the 
first and second term will cancel). In this regime the determination of Za from the 
RI/MOM condition is not feasible, unless one considers the p — )• oo case. As we have 
seen in sect. 14.21 in this limit these Goldstone pole terms become negligible. 

The simplest way to see a discrepancy between Ward identities and RI/MOM is 
through Zs/Zp ( jGiusti and Vladikas, 2000[ ). We start from the Ward identities ((X^ 
and p.34p . which we Fourier-transform to momentum space; cf. eqs. (j2.17p and (|2.28p . 
We next amputate the resulting correlation functions as in eq. (|2.29p and project them 
as in eqs. (|2.30p and (|2.3ip . obtaining the following vector and axial Ward identities: 



(mo2 - moi)rs(p;moi,mo2) = '^^^ ^ ^ {pl ^02) - S ^(p;moi) 



(m 



PCAC 
2 



PCAC\ 



12 



In the mass degenerate limit these become 



)rp(p;moi,mo2) = TrrTr 5 ^{p;moi)+S ^(p;mo2) 



, (4.24) 
. (4.25) 



Tsip) 



12 
1 

12 



-Tr 



-Tr 



dS ^ {p; nip) 
drrio 

{p; mo) 



(4.26) 
(4.27) 



The vector Ward identity and the definition (P^^ imply that Ts{p) = r„ (p). 

Once again we see from this and eq. (|2.22[) that if Zg is determined by the RI/MOM 
condition, then compatibility with this Ward identity requires that the mass renor- 
malization satisfy Z,„ = Zg^. 

We now turn to the axial Ward identity (|4.27p . The inverse quark propagator may 
be considered a function of either bare quark mass, toq or m^^^^. We differentiate 
the above axial Ward identity w.r.t. m^'^^'^ , obtaining 



Tp (p) 



,PCAC 



12 



-Tr 



dS'^p-, 



,PCAC\ 



Qj^PCAC 



drriQ 



iTr 
12 



dS ^ ip; mp) 
dmp 



Z^ 
Zp 



12 



-Tr 



dS ^ (p; mp) 
dmp 



(4.28) 
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In the last step, the substitution of the derivative dm.o/dm^''^^'^ by the ratio Zs/Zp 
is a consequence of eq. p.45|) . With the aid of the LSZ expression (|4.13p at zero 
momentum transfer (g^ = 0) we easily confirm that the second term of the l.h.s. 
of eq. (|4.28p diverges in the chiral limit, due to the presence of a Goldstone pole. This 
term is essential for the cancellation of a similar contribution of the first term on the 
l.h.s.. Again as p ^ oo these contributions become negligible. 

The ratio of the above vector and axial Ward identities (j4.26|) and ()4.28p gives the 
following result for ZpjZs'- 



Zp _ Tp 



Zs JTlPCAC 



(4.29) 



Tp 9mPCAC 



But this is not what the RI/MOM renormalization condition (j4.6p gives for the same 
ratio: 

Zp Pc 

-f = ^ ■ (4-30) 

The two determinations differ by a factor which, as we have already discussed, is 
necessary in order to cancel Goldstone pole divergences in Fp and which becomes 
negligible in the limit — oo; in practice this limit is realized as >> Aqcd- The 
absence of this factor from the RI/MOM determination is commonly referred to as 
Goldstone pole contamination. 

It is also instructive to identify the origin of the Goldstone pole in the non- 
perturbative part of the quark propagator: 



S~^(p;m) = i^j:i{p^;m;pi^)+mJ:2iP^;m;ti^)+J:3{p;m;n^) (4.31) 

where Sfc (with k = 1, • • • ,3) are form factors. The functional form of the first two 
terms is dictated by general symmetry arguments; the form factors Ei and E2 may 
be calculated in perturbation theory. The last term is a non-perturbative form fac- 
tor, known to several orders in the OPE; to 0(l/p^) it is given by ( [Politzer, 1976[ 
Pascual and de Rafael, 1982]): 



S3(p; m) = gl K {^^)- + O (p-4) , (4.32) 

where K is a mass-independent, gauge- dependent factor, in which logarithmic di- 
vergences have also been absorbed. Note the non-perturbative nature of this term: 
it is proportional to the chiral condensate {ipip) and vanishes in the large scale limit 
— )- 00. Upon inserting the quark propagator expressions (|4.3ip and ()4.32|) into Ward 
identities (|4:26| and (027]), we finOll 



''We have swapped, in eq. 114.2611 the dependence from the bare quark mass mo with that of the 
subtracted bare mass m. 
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rs{p;m) = E2b;m) + + O (p-^) (4.33) 

rpb;'") = ;;^PCAC ^2 (ap,a™) + ^ ^pcAc p2 + ^ (p^^) = 



Zp 



j:2{ap,am) + gl + O {p'^) (4.34) 



m p^ 



Thus, with the aid of an axial Ward identity, we see that, to 0{l/p^) in the OPE, the 
correlation function Tp has a pole in the quark mass, proportional to the chiral sym- 
metry breaking order parameter (ip'tp)- On the contrary, using a vector Ward identity 
we see that such a pole is absent, to 0{l/p^), from the scalar correlation function F^lfl. 

As stated repeatedly, this non-perturbative contribution to Tp vanishes like 
at large momenta. At finite momenta however, Tp is not regular in the chiral limit and 
p^m^^^^ >> Aqqj-, must be enforced. In practical simulations it is usually difficult 
to satisfy this requirement. A remedy consists in fitting Tp, at fixed momenta p^, 
by A{p'^) + B{p^)/m. Once the fit parameters A{p'^) and B{p^) are determined, the 
Goldstone pole contamination is removed from Zp by imposing the RI/MOM condition 
on the subtracted correlation function A{p'^) ( [Cudell et ai, 1999l|Cudell et ai, 2001[ ). 

As we have shown to 0{p~^), Ts does not suffer from Goldstone pole contami- 
nations and therefore the RI/MOM determination of Zs should be reliable even at 
relatively low renormalization scales. The same is true for the ratio Zp/Zs, computed 
from Ward identities, which are completely free of these problems. Thus, a reliable 
evaluation of Zp consists in multiplying the ratio of Zp/Zs, obtained from a Ward 
identity, by Zs in the RI/MOM scheme. The ratio Zp/Zs may be computed from 
hadronic Ward identities, cf. eq. (|3.44p . Alternatively, Ward identities in momentum 
space, based on the correlation functions Ts and Tp have also been used with satis- 



factory results ( Giusti and Vladikas, 2000 1 



Besides the single Goldstone pole, other, more complicated non-pcrturbative effects 
are of course also be present. Some of these could be revealed by higher orders in 
the OPE, which involve higher-dimension condensates ( [Pascual and de Rafael, 1982[ 



as p2 >> Ac 



Lavelle and Oleszczuk, 1992 1. In any case, all these contributions ought to disappear 

QCD- 



4.4 The RI/SMOM scheme 

In practice it is not easy to satisfy the bounds imposed by the renormalization window 
of eq. (|4.2p . At least one fine example where everything falls in place has been pro- 
vided in the literature ( [Becirevic et ai, 2004[ ). However, several cases exist in which 
the renormalization window turns out to be quite narrow: upon attempting to compute 
RI/MOM renormalization parameters at very high momenta, so as to reduce contami- 
nation due to the infrared chiral-symmetry breaking effects, one is faced with the prob- 
lem of discretization errors, which start flawing the data, even in an C'(a)-improved 

*An instructive exercise consists in plugging-in the above OPE expressions for Fg and Vp on 
the r.h.s. of eq. 114.291 1 in order to explicitly confirm that they combine to give the ratio Zp/Zs- 
Contributions proportional to the chiral condensate cancel out. 
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setup. The current precision of simulations is such that the infrared non-perturbatve 
effects are a significant source of error. Clearly it would be advantageous to renormal- 
ize the operators in a scheme which is free of this problem, suppressing these effects by 
choosing kinematics without channels of exceptional momenta. One such choice, which 
for quark bilinears is characterized by very simple and convenient kinematics is the re- 
cently proposed variant of the RI/MOM scheme, called RI/SMOM (SMOM stands for 
symmetric momentum subtraction scheme) . Although the new scheme has been intro- 
duced in the framework of domain wall fermions ( Aoki et ai, 2008[Sturm et ai, 2009 ) , 



its definition and the advantages derived from it are independent of the regularization 
details. 

Essentially the scheme consists in the choice of new kinematics for the vertex func- 
tion Go{pi,P2) of eq. (|2.28p . In the standard RI/MOM scheme the renormalization 
conditions for quark bilinear operators are imposed on Green functions with the oper- 
ator inserted between equal incoming and outgoing momenta, satisfying (in Euclidean 
space-time) 

Pi = P2 = ; q^Pi - P2 ^ . (4.35) 

In the above momentum configuration, the momentum inserted at the operator is 
therefore q = so that there is an exceptional channel, i.e. one in which the squared 
momentum transfer is much smaller than the typical large scale fP. These kine- 
matics define an asymmetric subtraction point. The renormalization procedure for 
RI/SMOM is very similar, but with the incoming and outgoing quarks having differ- 
ent momenta, pi and p2, satisfying 

pI ^ pI =q^ = p^ ; q = Pl~P2^0 . (4.36) 

There are now no exceptional channels. An example of such a symmetric momentum 
configuration on a lattice of linear extension L is 

2tt 27r 
p, = —(0,1,1,0) ; P2 = —(1,1,0,0) , (4.37) 

where p = 2TT/L{nx,ny,nz,nt). With this choice of a symmetric subtraction point, 
renormalized quantities such as [Go]ri [Ao]r and [EoJr depend only on a single scale 
p^ . The amputated Green function Aq is obtained once more as in eq. (|2.29p , whereas 
for the projected amputated vertex function Tq we introduce the following set of 
projectors: 

Ps = I ; Pp ^15 ; Pt = ^7/^71/ 

Pv = ^ ; Pa^ \q^.l5 i ■ (4.38) 

The RI/SMOM renormalization condition is that of eq. (|4.6p . imposed on the vertex 
functions [FoJr (implicitly defined though the above projectors), in the symmetric 
momentum configuration of eq. (|4.36p . We will use the shorthand notation {p^} for 
these kinematics. 

Comparing with eq. (|2.31|) . we see that the RI/SMOM definitions of Fs, Tp and 
Tt are identical to those of the standard RI/MOM scheme, while Fy and F^ have 
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been redefined through new projectors. The reason behind these changes is that we 
must ensure that the new scheme is consistent with Ward identities. If we wish to 
maintain condition (j4.9[) for the wave function renormahzation Z^, we must modify 
the definitions of Ty and Ta, as imphed by eq. (|4.38|) . These modifications must be 
accompanied by a redefinition of the quark mass renormahzation: instead of condition 
()2.22p . defined through (|2.20|) . wc must impose, in a symmetric momentum configura- 
tion, that 

hm |[r:„(p)]R - iTrL[A^]R75l| =1 , (4.39) 

where r^(p) is defined in eq. (|4.11[) and in eq. (|2.29p . Following a line of reasoning 
similar to that of sect. 14. 3[ it has been shown that this RI/SMOM scheme is consistent 
with Ward identities ( Sturm et ai, 2009] ). Clearly these choices are not unique. It 
is also possible to maintain the standard RI/MOM definitions for Ty and Ta but 
then the wave function renormahzation condition must be modified in order to ensure 
compatibility with Ward identities ( |Sturm et ai, 2009D . This is a different RI/SMOM 
scheme. 

Compared to the standard RI/MOM, the RI/SMOM scheme displays an improved 
infrared behaviour. We have explicitly shown in sects. 14.21 and 1 4.31 how the Goldstone 
pole contaminates the chiral limit of certain correlation functions as their quark ex- 
ternal momenta become exceptional (pi = P2) and the vertex operator docs not inject 
any momentum (q = 0). In the RI/SMOM scheme the kinematics are arranged so 
that this situation is avoided. In fact, using Weinberg' s theorem it has been shown 
that for the asymmetric subtraction point, chiral symmetry breaking effects vanish 
like 1/p^ for large external momenta p^. On the other hand, with bilinear opera- 
tor renormahzation constants defined at a symmetric subtraction point (with non- 
exceptional kinematics) , unwanted infrared effects are better behaved and vanish with 
larger asymptotic powers which are 0{l/p^) ( |Aoki et ai, 2008 ). Numerical evidence, 
in the framework of domain wall fermion discretization, provide strong support to 
these arguments ( |Aoki et al., 2008D . 



5 



Twisted mass QCD (tmQCD), 
renormalization and improvement 



In recent years, there has been a newcomer to the family of lattice actions, known 
as lattice QCD with a chirally twisted mass term (tmQCD for short). It consists 
in a modification of the mass term of the lattice Wilson fermion action. The new 
regularization has several advantages, but these come at a price. On the positive side, 
the tmQCD fermion matrix is free of the spurious zero modes which plague quenched 
simulations with Wilson fermions at small quark mass, as well as dynamical fermion 
algorithms. Moreover, the renormalization pattern of certain physical quantities is 
much simpler with tmQCD. Finally, the bare mass parameters of the theory may be 
tuned in such a way that improvement is "automatic" ; i.e. physical quantities such 
as masses and matrix elements are free of C(a)-effects, without having to intorduce 
Symanzik counter-terms. The price to pay is loss of symmetry: parity and time reversal 
are only recovered in the continuum limit. Moreover, flavour symmetry is also affected 
and this leads to loss of degeneracy between some hadrons. An important example 
concerns pions: the neutral pion mass differs from that of the two charged pions. 
Since this is due to discretization effects, degeneracy is restored in the continuum. 
A complete presentation of tmQCD is beyond our scope. The interested reader may 
consult existing review articles ( |Sint, 2007[rShindler, 20081 ). Here we will concentrate 
on the renormalization and improvement properties of the theory. 



5.1 Classical tmQCD 

It is instructive to begin by writing down the tmQCD classical Lagrangcan density 
in the continuum. For simplicity we will consider QCD with two degenerate flavours 
(called up and down quarks). The Dirac spinor in flavour space is then given by the 
doublet X = d) and the fermionic Lagrangean density is defined as 

>Ctm = X P + mo + iHqT^-f5 X ■ (5-1) 

Compared to the familiar standard QCD Lagrangean density, we now have an ad- 
ditional mass term which, being a pseudoscalar, is "twisted" in chiral space. Appar- 
ently, this is not classical QCD, as the extra twisted mass term breaks parity and 
SU{2) flavour symmetry. However, this is illusory. To see this, we redefine fermionic 
fields through chiral transformations (chiral rotations) in the third direction of flavour 
(isospin) space: 



Classical tmQCD 43 



exp |_i - 75 T J X 



X exp - 75 T 



X ^- x' 

We also redefine the two mass parameters through spurionic transformations 

mo — > ttiq = cos(q!) TTiQ + sin(a) iiq 
/ig ^ /i^ = cos(q;) /iq — sin(Q!) mo . 

Under these transformations, the Lagrangean density transforms as folfows: 



^ Am = X 



+ Too + ili'qT'^lb 



X 



(5.2) 



(5.3) 



(5.4) 



This means that the theory is form invariant under the ehiral rotations (|5.2p , provided 
they are accompanied by the mass spurionic transfrormations (|5.3p . In other words, 
these changes of field variables and mass definitions do not change the content of the 
theory. In fact tmQCD is a family of equivalent Lagrangean densities, connected by 
ehiral and spurionic transformations. A member of this family of equivalent theories 
is standard QCD, obtained through a specific transformation angle a, chosen so that 
M<3 ^ = 0. 

The above argument may be reformulated by defining an invariant mass Afjnv and 
a twist angle uj: 



ml + 



tan(a;) 



mo 



(5.5) 



The reason for this terminology will be clarified shortly. The Lagrangean density may 
be rewritten as 



X 



P 



Minv cxp (ia;r^75) 



X ; 



(5.6) 



i.e. it is a function of a single mass parameter Mjnv and a dimensionless angle u. The 
transformation (|5.4I) of the Lagrangean density may also be written as 



■'tm X 

with the invariant mass and twist angle defined as follows 



ML 



il^a 



tan(cLi') 



m,' 



(5.7) 



(5.8) 







Using eqs. (j5.3|) . it is easy to show that = Alinv', i.e. the invariant mass is indeed 
invariant under spurionic transformations. In physical terms, this is the quark mass, 
which in tmQCD is seen to be a combination of both the standard mass too and the 
twisted mass fiq. Morevoer, the new twist angle uj' can easily be expressed in terms of 
the old twist angle uj and the rotation angle a: 



tan(aj') ~ tan(aj — a) 



(5.9) 



Again, we see that tmQCD may be regarded as a family of equivalent theories, 
parametrized by an invariant mass Minv and a twist angle uj. Starting with a specific 
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tmQCD theory (i.e. a given value of defined through tan(a;) = fiq/mo), we obtain 
standard QCD by performing chiral rotations with a transformation angle a = w, 
which brings us to w' = (i.e. /i^ ~ 0). 

It is worth noting two values of the twist angle which arc of special interest to us. 
First, as already explained, for chiral rotations with a = uj, we obtain standard QCD 
(i.e. Lj' = 0, ji'q ~ and mp = Minv)- In this case, the quark mass is carried entirely 
by the standard mass mf,. The second case of interest is when the chiral rotations are 
a ~ uj — tt/2 (i.e. co' = 7r/2, fj,'^ = Mi^v and ttiq = 0). In this case, the quark mass 
is carried entirely by the twisted mass /i^. This is known as the maximally twisted 
theory. 

Since QCD and tmQCD at the classical level are equivalent theories, they should 
share the same symmetries. This means, for example, that loss of parity due to the 
presence of the twisted mass term in the tmQCD Lagrangean density is only apparent. 
It is indeed straightforward to confirm that Lagrangean density (|5.ip is invariant under 
the parity transformations : 

x = {x°,x.) x''' = (a;", -x) , 
Uo{x) Uoix'^) , 

Uk{x) ^ Ukix'^^k)^ , 

X{x) -> 70 exp[iw75r^] x(x'^) , 

X{x) X(a;^) exp[iujjr,T^] 70 . (5.10) 

This means that in the classical tmQCD framework, parity is still a symmetry, albeit 
with modified transformations of the fermion fields (cf. eq ()8.ip ). It should come as no 
surprise that the new parity transformations of x and x involve the very same chiral 
rotations that connect tmQCD to standard QCD. 

Analogous results may be obtained for the vector (isospin) symmetry. It is not 
really lost, as the twisted mass term of eq. (|5.ip may lead us to believe at first sight. 
Rather, it is simply transcribed as follows: 

Xix) cxp [ - i^75r^] cxp [«y t"] cxp [1^75^^] xi^) 

X{x) x{x) exp [i-75T-^] cxp [ - iy t"] exp [ - i-75T^] , (5.11) 

with 0"- (a = 1,2,3) the three rotation angles. We denote this symmetry group as 
SU{2)'y. Recall that from the beginning tmQCD has been formulated for a mass 
degenerate isospin doublet. 

Axial transformations are transcribed in the following form: 

X{x) cxp [ - i-75T^] exp [«yT''75] exp [i-l^T^] x{x) 

to 0^ UJ 

xix) x{x) exp [i-l5T^] exp [«yr°75] exp [ - ^-751-^] • (5.12) 

It is easy to verify that these transformations are a symmetry of the Lagrangean 
density (|5.6p . when Minv = 0. 
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The field rotations, relating standard QCD and tmQCD also relate composite field 
operators, defined in the two theories. Since tmQCD (with twist angle w) and standard 
QCD (with twist angle ut' = 0) are related by chiral rotations (|5.2[) with a = lj, the 
following relations hold between the vector and axial currents of the two theories: 

[^;]tmQCD = cos(c^) [V^hco ~ e^^" sin(c^) [^Iqcd a,6=l,2 , 
[A^^hmQCB = cos(c^) [AI]qcb - e'"' sin(L.) [V^hcB a,6 = l,2 , 

[4]tmQCD - [AlhcD . (5.13) 

Recall that a, b are flavour indices. Similarly, for the pscudoscalar density and the 
isospin singlet scalar density 5*^ = xx obtain 

[ntmQCD = [nqCD a,6==l,2 , 

[-P^JtmQCD = COs(tj) [P^]qcD " ^ sin(w) [S'%CD , 

[^°]tmQCD = cos(c^) [S%cu ^ 2i sin(w) [P%cu ■ (5.14) 

Finally, the familiar PCVC and PCAC Ward identities, have the following form in 
the tmQCD formulation: 

[^;]tmQCD = -2/.,e3afc [p^J.^QCD (5.15) 
[^^It-QCD = 2mo [P"]tmQCD + [^ItmQCD (5.16) 

5.2 Lattice tmQCD 

So far we have achieved precious little! The relation between classical QCD and 
tmQCD is simply a change of fcrmionic field variables, accompanied by a redefinition 
of the masses. A change of variables cannot bring about new Physics. Thus, classical 
tmQCD is simply an intricate way of writing down QCD. It is upon passing over to 
the field theoretic formulation that the equivalence between the two theories is less 
trivial and this has important consequences, both in renormalization and improvement 
properties of many Physical quantities. 

The proof that standard Wilson fermion QCD and tmQCD are equivalent Field 
Theories will only be shown schematically here. For a detailed demonstration, the 



reader is advised to consult the original reference ( [Frezzotti et ai, 2001[ ). The starting 
point is that the equivalence in question proceeds through linear relations between 
renormalized Green functions of the two theories. As a first step, we consider QCD 
and tmQCD, regularized on the lattice with Ginsparg- Wilson fermions. The exact form 
of the lattice actions is not important in this discussion; for example in both lattice 



actions the discretization of p may be the Neuberger operator (Neuberger, 1998 a 



Neuberger, 1998 & Neuberger, 1998c). What is important is that chiral symmetry is 
not broken, owing to the Ginsparg- Wilson relation obeyed by the regularized Dirac 
operator. This implies that the chiral rotations (|5.2p map the lattice tmQCD bare 
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action to the lattice QCD one, just like in the classical case. On the other hand, the 
variation of composite operators transform members of a given multiplet among each 
other, plus some extra terms. It has been shown that the presence of these extra 



terms does not affect the main line of reasoning (Frezzotti et ai, 2001 1. Thus they 
will be ignored in this discussion. The fact that chiral rotations transform both the 
tmQCD action and the composite fields into those of standard QCD implies that bare 
Green functions of the two theories are related. For example, the insertion of the scalar 
operator 5"" in a Green function gives rise to the following relation 

(•••^°---)§S'd = cosH(...5°...)S^cD + ^BinH (...p3...)t^n'^CD , (5.17) 

where the ellipses stand for other operators, which have the same form in QCD and 
tmQCD (e.g. V^, etc.) and the superscript GW indicates bare Green functions, 
regularized in a Ginsparg- Wilson lattice framework. The subscripts QCD and tmQCD 
denote the mass term regularization (standard or twisted) in which these Green func- 
tions arc defined. 

The next step is to pass from bare to rcnormalizcd Green functions. Since chiral 
symmetry is preserved by the Ginsparg- Wilson regularization, members of the same 
chiral multiplet, such as and , renormalize with the same renormalization factor 
Zga = Zp. In a mass independent renormalization scheme, this factor is the same both 
for standard QCD and tmQCD. This is easy to understand, since the two theories only 
differ in the way mass terms are introduced in the action. The bottom line is that when 
both sides of the last equation are multiplied by the same factor Z50 = Zp, we obtain a 
linear mapping between renormalized quantities, calculated in two frameworks (QCD 
and tmQCD): 

(•••['5'"]r---)qcD = COs(w)(---[S'"]r, •••)tmQCD + « sin(w) (• • • [P^]r • • • )tmQCD • 

(5.18) 

Thus, based on a specific lattice regularization which respects chiral symmetry, we 
have shown the equivalence of QCD and tmQCD beyond the classical level. So far 
the line of reasoning has been similar to the one of the classical theories: a change 
of variables and mass definitions in the actions and the path integrals maps QCD 
lattice (bare) expectation values into those of tmQCD. Multiplicative renormalization 
subsequently ensures that the same relations hold for the continuum (renormalized) 
Green functions. Just like in the classical case, we do not expect anything interesting 
out of these changes of variables. 

Things become non-trivial however, once we realize that the last expression, being 
true in the continuum, does not depend on the regularization details. This is a con- 
sequence of the principle of universality, which is a generally accepted assumption. In 
particular, universality implies that eq. (|5.18p is also true, up to discretization effects, 
for renormalized (continuum) Green functions, obtained from any other regularization; 
e.g. lattice Wilson fermions. Since chiral symmetry is broken by the Wilson term, the 
renormalization patterns of bare operators (which in the continuum belong to the same 
chiral multiplet) are now very different. This has important consequences. For exam- 
ple, as will be discussed below, the renormalization of with Wilson fermions is fairly 
comphcated, and this renders the direct computation of the l.h.s. of eq. (|5.18|) rather 
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cumbersome. On the other hand, the rcnormahzation of with Wilson fermions is 
much simpler (be it in standard QCD or in tmQCD). One may tune the mass param- 
eters of tmQCD so that the twist angle is ui = tt/2 and the first term in the r.h.s. 
of eq. (|5.18p vanishes. Rather than computing the l.h.s. using standard lattice QCD 
with Wilson fermions, one may then compute the r.h.s. with the tmQCD lattice reg- 
ularization of Wilson fermions, which has simpler renormalization properties. There 
are several interesting cases, besides the one sketched above, in which it is preferable, 
from the renormalization point of view, to work in a tmQCD framework. 

Having explained the main idea, we now present Wilson fermion tmQCD in some 
detail. The standard Wilson action p.2|) is modified by the addition of the twisted 
mass term: 



^ r ^ ^ — — ). — — ). -| 



The difference between the classical case and the present formulation is that the Wilson 
term breaks several symmetries. We saw, for example, that in the classical case parity 
is preserved in the modified form of cqs. (|5.10|) . This symmetry is now broken by 
the Wilson term of the action. It does however survive, when combined with flavour 
exchange 

X{x) ^ -i xix"") 70 , (5.20) 

where the gauge-field transformations are the usual onet0. Equivalently, instead of 
flavour exchange we can combine parity with a flip of the twisted mass sign, 

X{x) i 70 x(a;'^) , 
X{x) ~i x(a;^) 7o , 

liq ^ - l^l-q , (5.21) 

leaving the action invariant. Analogous considerations can be made for time- 
reversal. The consequence of this loss of symmetry is that matrix elements such as 
(0|V^"^(0)|7r°), which, due to parity, vanish in the standard Wilson theory, are non-zero 
in tmQCD. Of course, since we have argued in the previous section that tmQCD is 
a legitimate rcgularization of continuum QCD, such matrix elements will vanish in 
the continuum limit. They are lattice artefacts proportional to (some power of) the 
lattice spacing. Also note that often we study the asymptotic behaviour of correlation 
functions such as {P{x) P^{0)) at large time-separations, by introducing a complete set 
of states between the operators. In the standard QCD case these would be pseudoscalar 
states; in tmQCD loss of parity implies that more states (e.g. scalars) are also allowed. 

We now proceed to examine continuous symmetries. The twisted vector symmetry 
of eq. (|5.1ip is hard-broken by the Wilson term. Axial symmetry (cf. eq. (|5.12p ) is softly 
broken by the mass term Afinv as expected, but also by the Wilson term; the latter 

^The transformations of ea. lj5.20|l . with replaced by r^, are also a symmetry. 
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breaking is hard. Once the continuum and chiral limits are taken, these symmetries 
are expected to be restored. It is important to note however, that the subgroup of 
transformations 



obtained from cq. (|5.1ip by setting 9^ = 9^ ^ 0, remains a vector symmetry of the 
Wilson tmQCD action. We denote this group as U{\)y. Thus, due to the Wilson term 
we have the symmetry breaking pattern SU{2)y — > U(\)y. The reduced symmetry 
causes a lack of degeneracy between the neutral pion 7r° and the two degenerate charged 
pions TT^. The mass difference between charged and neutral pions is a discretization 
effect, proportional to (some power of) the lattice spacing, which vanishes in the 
continuum limit, where the full symmetry is restored. 

Another interesting symmetry is derived from the axial transformations (j5.12p 
which, upon setting a; = 7r/2 and 9^ = 9^ = 0, reduce to 



This group of transformations is denoted by U{1)\, though it has the appearance of 
the usual vector symmetry. Analogously, for 9^ = 9^ = we obtain a similar group of 
transformations, called U{1)'\. Indeed these are axial symmetries, being symmetries of 
the massless Wilson action, softly broken by the twisted mass term. Thus, in the chiral 
limit the twisted theory is symmetric under U{1)y (g) U{1)\ (g) U{1)\, which amounts 
to an SU{2) group, with one "vector" and two "axial" generators. It is not surprising 
that in the chiral limit, where the action (|5.19p reduces to the standard Wilson one, 
the full SU{2) symmetry is recovered. The interpretation of this symmetries as vector 
(in the standard case) or vector/axial (in tmQCD), depends on how the soft mass 
term is introduced in the action ( |Sint, 2007| ). 

5.3 Renormalization with tmQCD 

So far we have seen that the introduction of a twisted mass term in the Wilson fermion 
regularization breaks some discrete and continuous symmetries. We have also argued 
that this lattice theory goes over to QCD as the lattice cutoff is removed, so any 
effects due to loss of symmetry by the twisted lattice action vanish in the continuum 
limit. Having seen the shortcomings of tmQCD, it is high time we discuss some of 
the advantages. First of all, it is straightforward to see that the fermion determinant 
corresponding to the tmQCD action (j5.19p is positive definite, as long as /x^ 7^ 0. This 
is an important advantage for lattice simulations close to the chiral regime, but will 
not be further discussed here. A second advantage concerns renormalization. We have 
already argued, in rather general terms, that some operator renormalization is simpler 



Xix) exp [«yT"^] X{x) 
xix) x{x)exp [- jyr^] 



(5.22) 



X{x) exp [jyT-^] x{x) 




(5.23) 
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in tmQCD than in standard Wilson fcrmion lattice theory. Here we will discuss some 
important examples in detail. 

Before we do so, we should understand the renormalization properties of the two 
bare mass parameters of tmQCD, uiq and /z, (or equvalently, M-mv and w). The sym- 
metries of the tmQCD action suggest that, for degenerate flavours, the standard quark 
mass uiQ renormalizes as in eqs. p.l4p and (|2.15|) . while the twisted mass fiq renormal- 
izes multiplicatively. Moreover, the PCVC relation (|5.15|) is exact in lattice tmQCD, 
with the local vector current replaced by the point-split one p.lip . This in turn im- 
plies that the product of twisted mass Hq and pseudoscalar density (i.e. the r.h.s. 
of eq. (|5.15p ) is renormalization group invariant. Thus Zp^ renormalizes /igi 

[fiq]n = ^lq . (5.24) 

The above statements are generalizations of the arguments of Section |3l concerning 
standard lattice QCD with Wilson quarks. The reader should have no problem in 
convincing himself of their validity for tmQCD. 

Given these mass renormalizations, we define the twist angle through the following 
ratio of renormalized quantities: 

tan(a;) ^ Mi = ^ . (5.25) 

Note that with Ginsparg- Wilson fermions, chiral symmetry ensures that the two bare 
masses mo and fiq renormalize multiplicatively with the same renormalization factor, 
and thus it makes no difference whether the twist angle is defined through a ratio of 
bare or renormalized masses. It is the loss of chiral symmetry with Wilson fermions 
that leads to the above redefinition of the twist angle w.r.t. eq. (|5.5p . 

We stress that for mass independent renormalization, carried out in the chiral limit, 
there is no distinction between standard QCD and tmQCD, since the mass terms are 
absent from the action. Therefore, all renormalization constants computed in standard 
QCD simulations (e.g. Z,„, Zp) may also be used for tmQCD quantities. For instance, 
the renormalization factors of m and fiq in the last equation are the very same Zgo 
and Zp that renormalize the densities 5° and P in the standard Wilson theorjQ. 

With the above observations in mind, we can immediately generalize the tree level 
expressions (|5.13|) for the full theory. In particular, the continuum axial current A"^ 
(for a = 1,2) is expressed as linear combinations of renormalized tmQCD ones: 

[A^Jcont - Za [A;Jqcd = cos{uj) Za [A''^]tn.QCD + e'"' sm{uj) [V;']t,„QCD . 

(5.26) 

The subscripts QCD and tmQCD indicate bare quantities in the respective lattice 
regularizations. The first equation in the above expression links the continuum axial 
current to the lattice one in the standard formulation, as detailed in Section [3.21 The 
second equation does the same job in the tmQCD framework. As previously explained, 

^The same is true of the power subtraction rrici, being the counter-term of the standard bare mass 
mo- However, non-perturbatove determinations of ma, computed in a tmQCD framework, require 
extra care l |Frezzotti et al, 2006] . 
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the normalization factor Za is the same in both regularizations. The last term of the 
second equation contains the exactly conserved point-split vector current; thus the 
absence of a normalization factor. We may of course use, instead of [V^]tmQCD, the 
local current [V^]tmQCD, multiplied by Zy- Similarly for the pseudoscalar density 
(for a= 1,2) we have 

[-P"]cont = Zp [P^]qcU = Zp [P^tmQCD • (5.27) 

Equalities (j5.26p and (|5.27p arc valid up to discretization errors. 

Wc are finally ready to see a first advantage of maximally twisted tmQCD. To 
ensure maximal twist (i.e. uj = 7r/2), the standard bare mass uiq of the tmQCD 
regularization is tuned to its critical value nicr; cf. eq. (j5.25|) . Note that knowledge 
the ratio Zp/Zgo is not required for this tuning. Then, according to cq. ()5.26|) the 
axial current in the continuum limit is the vector current in tmQCD, which requires 
no normalization. Using eqs. (|5.26p and (|5.27p . we see that the standard PCAC Ward 
identity in the continuum corresponds to the PCVC relation 

K {V;:{x) P'(0))tn.QCD = (^'(^•) ^'(0))tmQCD • (5.28) 

X X 

This is an exact Ward identity in lattice tmQCD. Inserting, in standard fashion, a 
complete set of states between operators in the above expectations values, we obtain, 
in the limit of large time separations 

(0| [T>o(0)]tmQCD k) = 2^l, (0| [P(0)]t„,QCD k) . (5.29) 

Flavour indices have been dropped for simplicity. Combining these expressions and 
eq. ()5.26|) we get, for the pion decay constant 

U = — (0| [Ao(0)]co„t k> 

= hm — (0| [Vb(0)]tmQCD k> 

= hm ^ (0| [F(0)]tmQCD k> ■ (5.30) 

The last expression provides, within the tmQCD framework, a definition of which 
is free of any normalization or renormalization factors. Recall that in the standard 
Wilson fermion case, the computation of /^r from the matrix element (0|^ok) requires 
knowledge of the axial current normalization Za- Thus, the systematic error due to 
Za has been eliminated in tmQCD. 

Another example of the advantages of tmQCD is provided by the renormalization of 
the chiral condensate operator 5°. In the standard Wilson quark case, the symmetries 
of the theory imply the renormalization pattern (cf. eq. (|3.67p '): 

[5°]r = Zso \[S%cu + ^^] + ••• , (5.31) 



where the ellipsis stands for less vigorous power divergences (quadratic and linear) 
which, for dimensional reasons, are proportional to powers of the quark mass m (i.e. 
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they are absent in the chiral Umit). In tniQCD with maximal twist [u = 7r/2), this 
operator is mapped on (of. eq. (j5.14p ). The symmetries of lattice tmQCD with 
Wilson quarks imply the renormalization pattern 



[P^]t.„QCD + ^^^^ + ••• , (5.32) 



where the ellipsis again stands for less vigorous linear power divergences, which depend 
on the quark mass fiq. Thus, in the chiral limit the chiral condensate, regularized in 
tmQCD with maximal twist, is multiplicatively renormalizable. It is well known that 
power divergences are hard to control, whether calculated perturbatively or computed 
non-perturbatively (i.e. numerically). Suppose that the dimensionless coefficients cg 
and cp are calculated perturbatively. Such a calculation misses out dimensionless 
non-perturbative contributions, which are ©(oAqcd)- Since these contributions are 
to be multiplied by the power divergences and of the power subtractions, 
the pcrturbative calculation of cs and cp is inaccurate by terms which diverge in 
the continuum limit! On the other hand, if the coefficients cs and cp are calculated 
non-perturbatively, they are bound to be subject to discretization effects. These are 
expected to be 0(ajl. So again the error of the power subtractions diverges in the 
continuum limit. Nevertheless, it is remarkable that in tmQCD the chiral condensate 
suffers from less severe power divergences. If somehow one could work in the chiral 
limit of maximally tmQCD, determined non-perturbatively at very high accuracy (i.e. 
mcr ^ 0{a^)), then the power subtraction of eq. (|5.32p could be ignored altogether 
and the chiral condensate computation would be greatly simplified. 

These considerations have been extended to the calculation of matrix elements 
of dimension-6 four fcrmion operators, which determine the non-perturbative con- 
tributions in say, neutral Kaon oscillations and Kaon non-leptonic decays into two 
pions. In lattice QCD with standard Wilson quarks, the renormalization pattern of 
these operators is very complicated, due to the introduction of counter-terms which 
would have been absent, if chiral symmetry were preserved. In tmQCD it is pos- 
sible to map these operators into their partners of opposite parity, which have a 
much simpler renormalization pattern, in spite of the loss of chiral symmetry. Ex- 
plaining this in detail is beyond the scope of the present lectures. The interested 
reader is advised to consult the relevant literature, where these results are presented 
in considerable detail ( Frezzotti et al., 2001||Dimopoulos et al., 2006|Pena et ai, 2004 



[Frezzotti and Rossi, 2004 &[ )■ 



5.4 "Automatic" improvement in maximally twisted QCD 

Besides simplified renormalization patterns, tmQCD in its maximally twisted ver- 
sion has another nice property, commonly known as "automatic improvement" . By 
this we mean that physical quantities (e.g. masses, matrix elements etc.) have dis- 
cretization effects which are sub-leading; i.e. 0{a^). Improvement of these quanti- 
ties comes about by carefully tuning the twist angle uj = 7r/2, without the need 



^If Symanzik-improvement is used, these discretization effects are 0(a'^). 
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of introducing Symanzik counter-termfl Following the original proof of this state- 



ment (Frezzotti and Rossi, 2004a|, several simplified versions have appeared in the 
literature. In order to understand the line of reasoning of this section, the reader is 
strongly advised to acquaint himself with the relevant Sections of Chapter 2 of the 
present volume ( |Weisz, 2010[ ). 

The study of discretization effects of lattice observables is based on the Symanzik 
expansion, in the light of which the lattice action, close to the continuum, is described 
in terms of an effective theory 



S 



tm 



d-^y Co + a J d^y Ci + ■■■ . (5.33) 



As the notation in the above expression implies, our starting point is the tmQCD Wil- 
son fermion action (|5.19p at maximal twist (i.e. with mo = TOci)- Thus the Symanzik 
counter-terms reflect the symmetries of this action: 

>Co = x[]p + I K]r 75 ] X (5.34) 
^1 = -i [x cr ■ Fx] + [^^Jl [ XX ] ■ (5-35) 

Besides the above expansion for the action, Symanzik-improvemcnt implies that d- 
dimensional composite lattice fileds $iatt are also subject to an effective theory de- 
scription of the form 

$iatt = $0 + a $1 , (5.36) 

where $07 'I'l are continuum composite fields of dimensions d and d + 1 respectively. 
Actually, $iatt is a lattice transcription of $9. Combining eqs. (|5.33p and (|5.36p gives, 
the lowest-order Symanzik expansion for the vacuum expectation value of <I>iatt 

( $iatt)tm = ( $0 )o + a ( $1 )o - a J d^y {<Po Ci)o + ©(a^) . (5.37) 

The l.h.s. is the expectation value of the lattice operator in maximally twisted lattice 
QCD. The expectation values of the r.h.s. are continuum quantities; their subscript 
(•••)o indicates that thet are defined in terms of the continuum tree level action 
5*0 = J d'^yCo- Automatic improvement means that the last two terms on the r.h.s. 
vanish identically, due to some symmetries of the maximally twisted continuum theory. 



To prove this, we first explain what the relevant lattice symmetries are. We define 
the discrete chiral transformations in the Erst isospin direction TZl^ as: 



X ^- i 15 X 



X ^ * X 75 . (5.38) 
Next we define the operator dimensionality transformations V: 



*In the realistic case of QCD with non-degenerate quark masses, the non-pertubative determination 
of Symanzik counter-terms is a non-trivial and intricate task l |Bhattacharya et al., 2006[ |. 
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^5 ^ [/^9 -Mg] 




+ + + 


-5-* 


- - + 




- - + 


^ Ea; X /^9T^75 X 


+ - 



Table 5.1 "Parity" of the various terms of the fermion tmQCD action under discrete sym- 
metries. 

X{x) exp[^] x{~x) , 

X{x) x{-x) cxp[^] 
U^{x) -> Uli-x-afi) . (5.39) 

The gauge lattice action is invariant under the above transformations of the Unk fields 
Ufi{x). Composite operators of even (odd) dimension d arc even (odd) under these 
transformations, up to the flip of sign of the space-time argument x. Finally we define 
the twisted mass sign Rip transformation 

~ Hg (5.40) 

Each term of the lattice tmQCD action is even or odd under these transformations; i.e. 
it has definitive "parity" with respect to these symmetries. In Table \EA] we list these 
properties explicitly. It is then clear that the tmQCD fermion action is invariant under 
the combined transformations IZl^ (^"D® [/i^ — >■ — /Xg], which is therefore a symmetry 
of the lattice theory. On the other hand, the continuum tmQCD action Sq = J d'^yCo 
has positive 7?.5-parity (cf. eq. (|5.34l) and the first and last rows of Table [STT|) . Thus 
TZ^ is a symmetry of the continuum theory. 

Let us now examine each term of cq. (|5.37p in the light of these symmetries. It 
is important to keep in mind that the three vacuum expectation values on the r.h.s. 
are continuum quantities, determined by the continuum tmQCD action So = J d'^yCo, 
which is symmetric under TZ^. 

1. $iatt: Without loss of generality, we assume this operator to be even under TZl 
and to have even dimension d. Then its vacuum expectation value ( $iatt)tm is 
invariant under the combined transformations TZl^ ® T> [/Xg — ^ ~^^q\^ which is 
a symmetry of the lattice theory. This implies that the terms on the r.h.s. of 
eq. (j5.37p are also invariant under the same combined transformations. 

2. $0- This operator has positive T^-g-parity and even dimension d, being the contin- 
uum counterpart of the lattice operator $iatt- 
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3. $i: This operator has dimension (i+ 1, so it is odd under V. Since its expectation 
value ( 4>i)o is even under TZl ® V ® [fig ^ the operator i>i must be odd 
under TZl- This impHes that ( )o vanishes, being the expectation value of an 
7^5-odd observable $1, weighted by an T^-g-even action So- 

4. Ci: This is the C'(a)-counter-term of the continuum Lagrangean Cq- It is a con- 
tinuum dimension-5 operator. From eqs. (|5.35p and (|5.38p we deduce that Ci is 
odd under 7?,g transformations. 

5. $0 ^1- This is a product of two composite operators. The operator $0 is TJ-g-even 
while the operator Ci is T^-g-odd. It follows that the ( $0 ^1 )o vanishes, being 
the expectation value of an T^-g-odd observable $1, weighted by an T^g-cven action 

We have shown that the lattice expectation value ( <I'iatt)tm of eq. (|5.37p is equal to the 
leading order term of the Symanzik expansion ( $0)0 up to 0{a^) counter-terms. The 
0{a) counter-terms of the expansion are identically zero. Therefore, without introduc- 
ing Symanzik counter-terms, we find that expectation values of lattice operators in 
the maximally twisted theory are C'(a)-improved. Note that it is the vacuum expecta- 
tion values, rather than the maximally twisted action itself, which are automatically 
improved. 

The previous arguments have overlooked the following subtlety. The proof rests on 
the vanishing of the continuum vacuum expectation values ( $1 )o and ( $0 ^1 )oi 
due to their breaking of the discrete "chiral" symmetry TZ^, which is a symmetry 
of the continuum tmQCD action Sq- But do these "chiral condensates" vanish in a 
theory which, being a regularization of QCD, is characterized by spontaneous sym- 
metry breaking? The answer is affirmative, because the term generating spontaneous 
symmetry breaking is the twisted mass term, while possible 0{a) counter-tems are 
generated by the "chirally orthogonal" Wilson term, which in tmQCD breaks vector 
symmetry. In other words, the 7?.g-symmetry on which the automatic improvement is 
based is not really a chiral symmetry. It is a discrete subgroup of the flavour symme- 
try SU{2)v of eq. (|5.1ip . for w = tt/2 (i.e. maximal twist) and transformation angles 
(ai,a2,a3) — (0, 7r,0). Flavour symmetry is not spontaneously broken in the contin- 
uum, and thus the quantities ( $1 )o and ( $0 ^1 )o do indeed vanish by symmetry 
arguments. The situation is not altered by simulations, provided we take extra care 
that the continuum limit is approached before the chiral limit. In this way the chiral 
phase of the vacuum is driven by the mass term and not by the Wilson term. This is 
ensured by imposing the condition Hq » aAQ^^j^. 

We have shown that tmQCD is a variant of Wilson fermion regularization which, at 
the price of sacrificing some symmetries, enjoys simpler renormalization properties of 
composite operators and automatic improvement. The advantages on renormalization 
patterns are typically obtained by tuning the twist angle to a value which maps the 
original operator to its counterpart of opposite parity (e.g. axial current to vector 
current, scalar density to pseudoscalar density etc.). In some cases, this twist angle 
turns out to be different than 7r/2, the maximal twist value which ensures automatic 
improvement. Thus the renromalization and improvement advantages are not always 
satisfied simultaneously. Such difficulties arise with the four-fermion operators related 
to neutral Kaon oscillations and K tttt decays ( |Pena et al., 2004[ ). A way out of these 
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problems has been proposed ( [Frezzotti and Rossi, 2004&| , based on mixed actions with 
maximahy twisted valence quarks of the so-called Osterwalder-Seiler variety. Once 
again, the gains are accompanied by certain disadvantages, but these are issues beyond 
the scope of these lectures. 



6 

Conclusions 



The present lectures focus on three specific topics: (i) lattice Ward identities; (ii) the 
non perturbative RI/MOM renormalization scheme; (iii) renormalization and improve- 
ment of twisted mass QCD. From the study of these topics, significant theoretical and 
practical advantages may be obtained for the renormalization and improvement of 
lattice composite operators. Throughout we have used the Wilson rcgularization of 
the fcrmionic action. However, at least as far as the first two topics arc concerned, the 
issues discussed go beyond a specific rcgularization. 

Once chiral symmetry is lost by the rcgularization, Ward identities ensure its re- 
covery in the renormalized theory, as the UV cutoff is removed. In the Wilson fermion 
case chiral symmetry is broken by an irrelevant operator. The consequences are im- 
mediately seen in the power subtraction of the quark mass, the finite normalization of 
the axial current, the different renormalization constants of operators belonging to the 
same chiral multiplet etc. In other regularizations, such as domain wall fermions, the 
loss of chirality only appears at the "practical" level (i.e. when the fifth dimension of 
the domain wall is not infinite in simulations). The resulting loss of chiral symmetry 
is less manifest in such cases, but its effects arc analogous to those of Wilson fermions 
(albeit quantitatively less important). 

The RI/MOM scheme has been presented in considerable detail. Emphasis has 
been laid on the theoretical aspects of the scheme, the conceptual problems and their 
resolution, without presenting any results, which may be easily dug out in the lit- 
erature. Again, the language used was that of Wilson fermions, but it is quite clear 
that non-perturbative operator renormalization is nowadays a requisite for all lattice 
rcgularization schemes. 

Finally, the variant of Wilson fermions known as twisted mass QCD has been dis- 
cussed. Only selected aspects of this formulation have been exposed. The fact that 
the tmQCD fermion determinant is positive definite at non-zero twisted mass /ig has 
only been touched upon, although it is of great importance to simulations close to the 
chiral regime. Moreover, we have limited our discussion to tmQCD with two light de- 
generate quarks, although it is possible to generalize the theory so as to include heavy 
flavours ( [Frezzotti and Rossi, 2004c| . This subject has many faces, but the most im- 
portant properties of tmQCD are arguably those connected to renormalization and 
improvement. Significant simplifications of the standard Wilson formulation are ob- 
tained, which are however accompanied by loss of symmetry, recoverable only in the 
continuum limit. This is hardly surprising, being yet another example of the well 
known fact that most significant theoretical gains in Quantum Field Theory come at 
a price. 
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Appendix: Spurionic chiral symmetry 



Fermion fields are decomposed into left- and right-components iJj ~ 4'l + ipRi with 

/ 1- 75 , 7 7 1 +75 

Vl = — ^ — ip , Wl = V — ^ — , 

/ 1 + 75 , 7 7 1 - 75 ,„ ,^ 

WR = — ^ — W , ipR = Tp — - — . (7.1) 

They transform under the chiral group SU{Nf)l ® SU{Nf)r as follows: 

V'fl^V'H ^ UriI^r , V^i?, ^Vifl = V^i?, t/j, , (7.2) 



with 

f^L,H = exp 



l,...,iV2-l (7.3) 



l,r y 

The vector transformations ()2.6p are recovered for C/l = C^_r (i.e. = a^), while the 
axial transformations (j2.8p are recovered for C/^ = J/jj, (i.e. af^ = —ol^)- 

The kinetic term of the lattice action (|2.2p is invariant under these transformations, 
while the mass- and Wilson-terms break chiral symmetry. However, the symmetry may 
be enforced on the mass term, once the mass matrix Mq is generalized, so that it is 
neither real nor diagonal (the physically relevant case corresponds to Mg — Mq and 
diagonal) . In order for the action to remain Hermitean, the mass term must be modified 
as follows: 

tpLMoi^R + i^RMoi^L ^ i^LMoiJjR + ^rMIiPl ■ (7.4) 
We also impose the following chiral transformation of the mass matrix: 

A/o ^ - Ul Mo . (7.5) 

Transformations of the mass parameters arc called spurionic. Under the transforma- 
tions (|7.2p and (|7.5p the modified mass term of the action is chirally invariant. Note 
that the Wilson term in the action remains a (spurionic) chiral symmetry-breaking 
term. 

The generalized, non-diagonal bare mass matrix Mq may be decomposed into 
non-singlet and singlet quark mass contributions (this terminology will become clear 
shortly): 

Nl-l 

Mo = ^ to'' A" + m'^^/ . (7.6) 
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Using the trace property of eq. p.7p for the group generators A" we see that 

= ^Tr[A/oA''] , = ]^ Tr[A/o] • (7.7) 

Clearly m'^^ is the average of the diagonal elements of Mq . It is straightforward to show 
that under infinitesimal chiral transformations these mass components transform as 
follows: 

. r_ avl ' _aV| ^ / a a \ ~a o\ 

m — Im I = m + (a^ — ol^) m . ( ' -oj 

In the special case of vector transformations {a1 = af, = a°), the above simplify to: 

to'^v ^ ^^avj' ^ ^av ^ (7 g^ 

Thus, the non-singlet mass components are indeed SU{Nf)v multiplets in the 
adjoint representation (triplets for Np = 2, octets for Np = 3 etc.). On the other hand, 
m'^^ is a flavour singlet. Under infinitesimal axial transformations {a1 = —a% = a"), 
eqs. (|7.8p reduce to 

m" [to'']' = m'^ + d"^" a^m" + ia" m"^ , 

TO^^ [to^^I' = TO^^ + — a" . (7.10) 

Therefore, flavour non-singlets m° and the singlet mF' transform into each other under 
axial transformations; i.e. they belong to the same chiral multiplet. 
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Appendix: Lattice discrete 
symmetries 



In this Appendix we gather the discrete symmetry transformations of the Wilson 
action , (j2.2|) ( [Bernard, 1989 1 . Although of rather limited value to the present set 



of lecture notes, they are generally useful for the classification of physical states, the 
mixing of operators under renormalization etc. We start with parity, which is defined 
as follows: 



Time reversal is given by 



(x°,x) 


-> 


x'^ = (a;°,-x) 


Uoix) 




Uoix'^) , 






Ukix'P-k)^ , 


tp{x) 


-> 






—?■ 


70 ■ 


(x°,x) 


-> 


x'^ = (-a;°,x) 






U„{x'^-6y , 


Uk{x) 




Ukix'^) , 


ip{x) 




70 75 ipix"^) , 






i^{x'^) 75 70 • 



Charge conjugation is defined as 



(8.1) 



■2) 



u^{x) ^ u^{xy , 

ip{x) C ij^{x) , 

'4>{x) -V{x) C-^ . (8.3) 

The T superscript denotes transpose vectors and matrices. The charge conjugation 
matrix satisfies 

C 7m C - -7; - -7J • (8.4) 
A realization of C is C = 7072 • 
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Appendix: Regularization dependent 
scheme 



The acronym "RI" stands for "regularization independent" and the reason for this 
is mostly historical: before the advent of lattice non-perturbative renormalization, 
Zo-factors were calculated in lattice perturbation theory, with MS as the scheme of 
preference. This means that two perturbative calculations had to be carried out; one 
in the continuum and one on the lattice. Schematically, in dimensional regularization 
(DR) the continuum perturbative expression for the bare correlation function Tq, at 
one-loop is 



1 



(47r)2 



(0) 



1 



(0) 

Tr 



1 (P' 



) + Cr™ ) 



(9.1) 



where 7p'''' is the LO anomalous dimension of the correlation function Tq- To this 
order it is a universal, scheme-independent quantity. As implied by the notation, the 
finite constant Cj? ^ depends on the regularization scheme and the chosen gauge. The 
factor 1/e is an abbreviation for 

\ = I + ln(47r) - 7E , (9.2) 

where e = {4 — D)/2 and 7^; stands for Euler's constant. In DR we work in D di- 
mensions, where the original bare coupling go has dimension e. Here the scale fi is 
introduced to render the bare coupling 50 (m) dimensionless. The /^-dependence of the 
r.h.s. of eq. (j9.1|) is only apparent. 

Imposing the MS renormalization condition amounts to removing the 1/e diver- 
gence. Since the renormalization constant of the projected Green function Tq is given 
by Zr = Z^^Zq (see eq. (|2.36p ). this implies for Zr the value 



7MS,DR 



(5o(/^),e) = 1 - 



.90 (/^) 
(47r)2 



(0) 1 
Tr 7 



(9.3) 



Consequently, the renormalized Green function is given by 



lim 

e-i-O 



= 1 



(47r)2 



(9.4) 



where, to this order, we are free to replace go by ff5jg(/t), the MS renormalized coupling 
constant. 
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The same calculation can be repeated on the lattice. Now the UV cutoff is provided 
by the inverse finite lattice spacing and thus the 1-loop calculation yields: 

r^^^b,go(a),«) = 1 + ^ ( - 7^°^ Hpaf + Cp^^^) + 0{a) , (9.5) 

where go (a) is the bare coupling of the lattice action. The renormalization scheme 
can again be chosen at will; as previously stated, the MS is often chosen also on the 
lattice^. Thus, we must satisfy the renormalization condition 



lim 

a-)-0 



(9.6) 



where again the lattice coupling .go(i) should be traded for the MS renormalized 
coupling constant gy[g{fJ-) ■ This point of principle is of limited relevance for a 1-loop 
calculation. From eqs. (|9.4p . (|9.5p and (|9.6p . the following renormalization constant is 
obtained: 



Zp {fia,go{a)) 



1 



.90 (q)^ 
(47r)2 



-fLAT 



(9.7) 



We now recall that the renormalization constant of the amputated vertex To is 
Zr = Z^^Zo- The quark field renormalization Z^ can be calculated, from the quark 
propagator S{p), with an analogous procedure; c.f. eq. p.22p . The result is 



,MS,LAT 



{fia,go{a)) = 1 + 



.90 (a) ^ 



Combining eqs. (|9.7p and (|9.8p we obtain 

50 (a) ^ 



7^°^ln(Ma)^ + a 



2 I /-(DH _ C^^^ 



,MS,LAT 



{fia,go{a)) 



1 



(47r)2 



7g"ln(^a)2 + Ap + As 



(9.8) 



(9.9) 



where 



7o = 7r + 7s , 

Ap = - Cj^AT , (9.10) 
s 



A r^DR ^LAT 

As = Cs - C- 



It is this renormalization constant (with this choice of renormalization condition) which 
is usually denoted by Zq in lattice perturbation theory calculations. The dependence 



^This seemingly unnatural choiee (the MS is closely linked to continuum DR) has a few advantages. 
For example, matrix elements of effective Hamiltonians, once calculated non-perturbativcly on the 
lattice, must be renormalized and combined with perturbatively calculated Wilson coefficients, in 
order to obtain physical amplitudes; cf. eqs. I|4.3|l and l|4.4|l . The renormalization-group invariance of 
these amplitudes is guaranteed only if the Wilson coefficients and the renormalization constants are 
calculated in the same renormalization scheme. Since the former are often known in the MS scheme, 
this scheme is also preferred for the calculation of the latter. 
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of Z, 



MSXAT 



on the coefficients Cj? ^ and Cl?^ comes from the choice of the MS renor- 



mahzation condition (see eqs. (|9.4p and (|9.6p ) whereas its dependence on Cp^""^ and 
C^^""" from the lattice regularization (see eq. (|9.5p ). Two perturbative calculations are 
thus necessary, one in the continuum for the C'^^'s and one on the lattice for C^^'^^s. 
The presence of the C'^^'s on the r.h.s. of eq. ()9.7p is sometimes referred to as the 
"regularization dependence" of the renormalization scheme. 

On the other hand, the RI/MOM scheme is obtained by imposing condition (|4.6p 
on the lattice correlation function (|9.5p . It follows that the renormalization constant 
of interest, in one- loop perturbation theory, is given by 



, RI/MOM 



(//a, 50(a)) = 1 + 



(47r)2 



Similarly, condition (j2.22p . related to the quark field renormalization, gives 



(9.11) 



7RI/M0M 



(A^a, 30(a)) 



1 



3o(a)^ 



The last two expressions combine to give 



a 



LAT 



(9.12) 



„Ri/MOM. , , , 5o(a)^ 

Zq (Aia, 50(a)) = 1 + j^;^ 



-.LAT 



, (9.13) 



Clearly, the above RI/MOM result does not depend on a continuum regularization 
(such as DR), but only on the lattice. This is referred to as "regularization indepen- 
dence" . 
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